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Abstract 


The  directions  and  the  differences  of  the  principal  membrane  and  bend¬ 
ing  stresses  at  a  point  of  a  shell  or  plate  are  calculated  from  the  bire¬ 
fringence  observed  at  normal  incidence  on  two  birefringent  coatings,  one  on 
each  face  of  the  shell .  A  unique  neutral  surface  is  assumed ,  which  is  exact 
for  identical  coatings  or  for  identical  Poisson's  ratios  in  shell  and  coat¬ 
ing.  Explicit  formulas  are  obtained  with  the  simplifying  approximations  of 
a  negligible  effect  of  the  rotation  of  the  principal  directions  and  of  a 
linear  variation  of  thq  magnitude  of  the  principal  stress-difference  over 


the  finite  coating  thickness.  These  assumptions  are  strictly  valid  for  very 


thin  coatings,  but  give  also  reasonably  good  results  for  ordinary  birefrin¬ 


gent  coatings. 
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The  validity  of  the  derived  expressions  was  checked  experimentally  in 
square  plates  subjected  to  anticlastic  bending  and  membrane  stress  of  known 
magnitude  at  various  orientations.  Two  types  of  tests  were  performed.  At 
first  [1]*  the  membrane  tension  was  applied  through  a  series  of  pins  in  holes 
around  the  perimeter  of  a  plate  with  relatively  thick  and  optically  insensi¬ 
tive  coatings.  Agreement  between  observed  and  calculated  birefringence  under 
these  unfavorable  conditions  was  reasonably  good  at  the  higher  loads,  but  not 
good  enough  to  permit  an  accurate  inverse  calculation  of  the  stress-differences 
from  the  observed  birefringence.  The  errors  were  found  to  be  partly  due  to 
the  difficulty  of  applying  a  true  membrane  force  through  pins  without  intro¬ 
ducing  any  bending. 

To  assess  accurately  the  new  method,  the  troublesome  pin  loading  method 
was  abandoned  and  the  membrane  stress  was  simulated  by  a  residual  stress  dis¬ 
tribution  in  the  coatings.  The  initial  method  of  applying  anticlastic  bend¬ 
ing  was  retained.  Three  plates  were  tested,  a  trial  one  with  0.164  in.  thick 
coatings  and  two  more  with  0.108  *  0.0015  and  0.057  *  0.0015  in.  thick  coat¬ 
ings.  At  first  the  observed  birefringence  was  compared  with  corresponding 
values  calculated  from  the  applied  bending  and  membrane  stresses.  The  real 
test  of  the  proposed  new  method,  however,  was  the  inverse  calculation  of  the 
principal  bending  and  membrane  stress  differences  from  the  birefringence  ob¬ 
served  on  the  two  coatings.  The  agreement  with  the  applied  values  was  quite 
good,  especially  at  the  more  interesting  higher  loads,  even  when  a  large  ro¬ 
tation  of  principal  stress  occurred  within  one  of  the  coatings. 


List  of  Symbols 


A,  I 


Effective  cross-sectional  area  and  moment  of  inertia 
per  tinit  width  of  composite  shell 


C 


(T 


Stress  optical  coefficient  in  lb/(in  x  fringe) 


E*,  E 


Young's  moduli  of  shell  and  coating,  respectively 


M,  »  M, 


Major  and  minor  principal  bending  moments  per  unit 
width  of  shell,  positive  when  causing  tension  at  z  >  0 


M  =  3!  0  Principal  difference  of  bending  moments 

N  i  N-  Major  and  minor  principal  membrane  forces  per  unit 

width  of  shell 


N  =  -  Nj  ^  0  Principal  differences  of  membrane  forces 

R  =  n/8  Ratio  of  the  rate  of  rotation  of  the  principal  stresses 

in  the  wave  front  to  half  the  rate  of  change  of  the  rela¬ 
tive  phase  difference  at  a  point  of  the  light  path 


real  rotation  through  the  coating  to  half 
the  retardation  for  single  passage  calculated  from  the 
real  stress  at  mid-thickness. 


(1+r2) 


1/2 


a.  >  0 
1 


Distance  from  centroid  of  composite  section  to  mid¬ 
thickness  of  coating  above  (i  =  o)  or  below  (i  =  u) 


2H  >  0;  2h  >  0 


Thickness  of  coated  and  uncoated  shell,  respectively 


h.  >  0;  H.  >  0 
1  *  1 


Ah.  >  0 
1 

i  0 


m 


n^£  i  Oj  i  0 


bi 


n 


in  . ;  n  • 

;  *  mi  ’  * 


all  i  0 


z 


Distance  from  centroid  of  composite  section  to  inter¬ 
face  and  to  free  surfaces  of  coatings  respectively, 
above  (i  =  o)  or  below  (i  =  u),  (Fig.  1) 

Thickness  of  coating  above  (i  =  o)  or  below  (i  =  u) 

o 

Fringe  order  (X  =  5461  A)  observed  in  reflected  light 
on  coating  above  (i  =  o)  or  below  (i  =  u) 

(n^  -  n^)/(n^  +  n^) 
o  u  o  u 

Fringe  order  observed  above  (i  =  o)  or  below  (i  =  u) 
when  shell  is  subjected  only  to  membrane  forces  N^,  N^, 
or  only  to  bending  moments  M^,  M^,  respectively 

Fringe  order  which  would  appear  under  both  M  and  N  if 
the  stress  throughout  the  coating  were  equal  to  the  total 
stress  at  a  distance  z ;  at  the  coating  mid-thickness  above 
(i  =  o)  or  below  (i  =  u);  or  at  a  distance  z^ 

Distance  from  centroidal  surface  of  composite  shell,  posi¬ 
tive  downwards 


a  Angle  from  the  major  principal  membrane  force  to 

the  major  bending  plane  (M.),  ccw  (counterclockwise) 
when  viewed  from  below  (z  >  0) 

2B  Rate  of  change  of  relative  phase  difference  over  a  unit 

length  due  to  birefringence  with  no  rotation 

Angle  from  the  reference  direction  to  the  major  princi¬ 
pal  membrane  force  and  to  the  major  bending  plane 

measured  ccw  positive  when  viewed  from  below  (z>  0) 

6  %  0  MAH/NI  =  (n,  ./n. . )H/a.  i  0  (same  for  i  =  o  or  i  =  u) 

bi  ti  1 

O.j,  0^  Angles  from  the  reference  direction  to  the  major  princi¬ 

pal  total  stress  at  the  surface  of  coating  above 
(i  =  o)  or  below  (i  =  u),  or  at  distance  z  from  the  neu¬ 
tral  surface,  respectively.  All  angles  measured  ccw 
positive  when  viewed  from  below  (z  >  0) 

Poisson's  ratio  of  shell  and  coating,  respectively 

Major  and  minor  principal .bending  stresses  in  the  coating 
at  distance  z 

%z  ~  *^lz  ”  *^2z  ^  0, Difference  of  principal  bending  stresses  at  a  distance  z 

and  at  mid-thickness  z  =  a.  of  coatings  above  (i  =  o)  or 
^  ^  Q  below  (i  =  u).  Their  major  principal  directions  are  par- 

bi  '  allel  to  the  corresponding  major  principal  stress,  hence 

change  by  90®  during  passage  from  lower  to  upper  coating 

Major  and  minor  principal  membrane  stresses  in  the  coating 

Og  i  0  Difference  of  principal  membrane  stress.  Its  major  prin¬ 
cipal  direction  is  parallel  to 

Major  and  minor  total  principal  stress  at  distance  z 

Difference  of  principal  bending  stress  at  a  distance  z. 

Its  major  principal  direction  is  parallel  to 

Angle  from  the  major  total  stress  a.  to  the  major  bending 
]b 

stress  at  the  plate -coating  interface  (1)  and  the  free 
surface  (2)  of  coating,  above  (i  =  o)  and  below  (i  =  u), 
or  at  distance  z,  respectively.  Angles  are  counterclock¬ 
wise  positive  when  viewed  from  below  (z  >  0) 

n  Rate  of  rotation  per  unit  length  of  the  principal  stresses 

in  the  wave  front,  positive  when  left-handed 


"z  =  “iz  -  °2z  »  “ 
*il’  *12  ’  *z 
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Introduction 

The  determination  of  the  surface  strains  of  opaque  bodies  by  the  use  of 
a  cemented  birefringent  coating  and  reflected  polarized  light  at  normal  in¬ 
cidence  has  been  suggested  by  Mesnager  [2]  since  1930  and  studied  by  Oppel  C3] 
(1937).  The  method  was  first  developed  to  practical  use  at  Brown  University 
[4-5]  over  thirteen  years  ago,  and  independently  in  France  [6-7]  and  later  in 
Japan  [8].  This  was  made  possible  by  the  use  of  the  highly  strain-optically 
sensitive  epoxy  resins  [5],  which  firmly  adhere  to  many  solids  and  may  show 
several  fringes  even  at  a  relatively  small  thickness  and  at  strains  as  snail 
as  those  of  metal  structures.  The  method  allows  the  determination  of  the 
maximum  in-plane  shear  strain  at  all  points  of  the  surface  from  a  single 
picture  almost  as  easily  as  in  ordinary  photoelasticity.  This  is  exact  for 
regions  with  small  strain  variations  and  far  from  the  coating  boundaries. 

In  the  presence  of  strong  strain  gradients  and  of  curving  of  the  sxirface, 
the  strains  vary  through  the  thickness  of  the  coating.  Considerable  errors 
[9]  may  result  if  the  analysis  is  based  on  uniform  strains  over  the  thick¬ 
ness  of  the  coating.  In  plates  and  smooth  shells,  strong  variations  of  cur¬ 
vature  or  of  strain  in  the  metal-plastic  interface  over  short  distances  (of 
a  few  coating  thicknesses)  are  unlikely  and  the  above  errors  should  be  neg¬ 
ligible.  However,  another  difficulty  may  arise  from  the  superposition  of 
bending  and  membrane  stresses  in  the  coating  and  their  effect  on  polarized 
light.  Pure  bending  produces  a  linear  strain  variation  through  the  coating 
thickness  which  can  be  easily  calculated  from  the  resulting  birefringence 
[14].  Pure  membrane  strains  can  also  be  directly  determined  from  their  op¬ 
tical  effect  since  they  are  constant  over  the  coating  thickness.  But  the 
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superposition  of  bending  and  membrane  loading  at  different  principal  direc¬ 
tions  produces  principal  stresses  which  vary  in  intensity  and  direction 
across  the  coating  thickness.  The  photoelastic  effect  in  such  variable 
stress  fields  is  very  complicated  ClO-25]  and  no  exact  inverse  solution 
for  the  stresses  in  terms  of  the  total  optical  effect  has  been  obtained. 

Coatings  of  infinitesimal  thickness,  however,  would  obviously  show  no 
rotational  effects,  and  very  thin  coatings  would  probably  show  very  little. 

It  appears  interesting  therefore  to  inquire  whether  in  coatings  of  the 
usual  thickness  (0.050  to  0.250  in.)  the  total  relative  retardation  may  not 
also  be,  to  an  acceptable  approximation,  proportional  to  the  integrated 
principal  stress  difference  independently  of  rotation.  The  determination 
of  bending  and  membrane  stress  from  the  total  birefringence  would  then  be 
relatively  simple,  and  hopefully,  sufficiently  accurate.  Obviously  some 
modifications  are  needed  for  the  determination  of  the  directions  and  differ¬ 
ences  of  both  principal  membrane  forces  and  of  principal  bending  moments  at 
a  point  of  a  plate  or  shell.  This  is  the  only  information  which  can  be  ob¬ 
tained  from  normal  incidence  measurements,  which  are  independent  of  any  in¬ 
plane  isotropic  state  of  stress.  The  individual  principal  components  may  be 
afterward  determined  with  the  help  of  additional  observations  at  oblique  in¬ 
cidence,  or  by  some  interferometric  measurement  of  absolute  retardations  [33], 
or  ,  for  plates  and  thin  shells,  from  the  principal  stress  differences  and 
directions  by  a  method  using  integration  from  a  boundary  separatley  for  mem¬ 
brane  and  for  bending  stresses,  as  suggested  by  Akhmetzyanow  [26]. 

Accordingly  the  main  problem  is  to  find  the  principal  membrane  and  bend¬ 
ing  stress  differences  and  directions.  Four  items  of  information  are  re¬ 
quired  to  determine  these  four  unknowns,  and  can  be  obtained  from  normal 
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incidence  measvireinents  with  two  birefringent  coatings,  i.e,  from  two  relative 
retardations  and  two  principal  directions,  as  was  shown  earlier  [1]. 

A  similar  problem,  that  of  a  shell  with  a  reflective  sheet  at  mid- thick¬ 
ness,  has  been  recently  treated  by  Kayser  C23D  and  Kuske  [24].  They  consider 
the  rotation  of  the  stress  and  calculate  the  optical  effect  from  the  known 
stresses  for  a  large  number  of  cases.  The  inverse  problem  is  then  solved  by 
comparison  with  the  obtained  direct  solutions.  The  method  would  not  be  con¬ 
venient  for  birefringent  coatings  as  many  direct  solutions  should  be  obtained 
for  every  ratio  of  metal  to  coating  thickness.  An  approximate  inverse  solu¬ 
tion  for  shells  with  thin  birefringent  coatings  (coating  about  one  tenth  of 
plate  thickness)  has  been  recently  made  by  Aben  [25],  on  the  basis  of  the 
differential  equations  of  birefringence  in  inhomogeneous  stress  fields  from 
observations  with  light  of  different  wavelengths. 

Analysis  of  the  Stresses 

Only  the  more  likely  method  using  coatings  on  both  faces  will  be  studied. 
The  thickness  of  the  metal  shell  or  plate  is  2h  ,  of  the  coated  shell  2H 
and  of  the  two  coatings,  above  and  below,  Ah^  and  Ah^  respectively  (Fig.  1). 
The  rectangular  coordinate  system  has  the  axes  x  and  y  in  the  centroidal 
plane  and  the  z-axis  perpendicular  to  it  and  pointing  downward  (z  >  0  below). 
The  conditions  of  continuity  of  the  tangential  strains  across  the  interface 
lead  to  the  conclusion  that  normal  in-plane  forces  per  unit  width  ^  on 
two  perpendicular  planes  give  also  rise  to  bending  unless  they  are  applied  at 
specific  distances  from  the  mid-surface  of  the  shell.  In  general  these  dis¬ 
tances  will  be  different  in  the  two  directions  and  will  depend  on  the  ratio 
Ni  ;  N2  ,  so  that  no  effective  centroid  exists.  Likewise  bending  moments  per 
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umit  width  )  M2  in  general  result  in  bending  about  two  different  neutral 
axes  in  the  two  directions  [28].  However  for  ^  N2  and  ^  ^2  ® 

tinct  neutral  surface  still  exists  when  the  two  coatings  are  of  identical 
thickness,  or  when  shell  and  coating  have  equal  Poisson's  ratios.  Then  the 
principal  membrane  and  bending  stress-differences  0^  and  at  a  dis¬ 

tance  z  in  the  coating,  defined  as  positive  quantities  with  principal  di¬ 
rection  parallel  to  the  corresponding  major  stress,  and  the  resulting  posi- 

^  lb 

tive  birefringence  n^  or  n^  respectively,  can  be  easily  found  from  the 
corresponding  positive  principal  differences  N  of  membrane  forces  and  M 
of  bending  moments  (Fig.  1). 

N  =  Ni  -  N2  i  0  M  =  -  M2  i  0  (1) 

=  N/A  5  0  -  |Mz/l|  =  aj^,lz/a.|  i  0  (2) 

where  is  the  bending  stress-difference  at  mid-thickness  of  coating  i  . 

Wherever  possible  formulas  are  given  in  terms  of  the  fringe  orders  n^  and 
which  would  be  observed  under  ptH>e  membrane  or  bending  loads,  but  can  al¬ 
ways  be  written  in  terms  of  N  and  M  with  the  following  substitutions 


n. .  =  20.Ah,/C  =  2NAh./AC  ^  0 

ti  t  i  0  10 


n,  .  =  20,  Ah./C  =  2Ma.Ah./IC  ^  0 
bi  bz  10  110 


(3a,b) 


where  the  subscript  i  may  be  either  o  (denoting  the  coating  above)  or  u 
(below),  but  the  same  in  any  one  equation.  The  denominators  A  and  I  rep¬ 
resent  the  effective  area  and  moment  of  inertia  per  unit  width  of  the  coated 
shell,  given  in  equations  (4a,  b)  below  in  terms  of  the  positive  distances 

h  a  n  d  h  from  the  neutral  surface  to  the  interfaces  and  a  ,  a  to 
o  u  o  u 
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the  raid-thickness  of  the  coatings,  the  positive  thicknesses  Ah^  and 

Ah^  o  f  the  corresponding  coatings  (Fig.  1)  and  the  Young’s  moduli  and 

Poisson's  ratios  E*  ,  E  and  v*  ,  v  of  shell  and  coatings  respectively. 

Expressions  (4a,  b)  are  exact  when  either  Ab  -  Ah  or  v*  =  v,  or  when 

o  u 

=  N2  and  .  The  siibscripts  o  (for  over)  and  u  (for  under) 

indicate  the  side  from  the  centroidal  surface. 


* '  iri-w  ‘”0  *  ''u>  ♦  '*0  *  **'u 


\  (Ha,b) 


'  =  *  il  ("-o’'’ " 


In  general  v*  /  v ,  Ah^  j*  Ah^  ,  neutral 

surface  exists.  These  values  are  then  incorrect,  but  may  be  considered  as 
good  approximations  because  the  differences  between  v  and  v*  ,  and  be¬ 
tween  Ah  and  Ah  are  usually  small.* 

The  total  principal  stress-differences  and  directions  at  a  distance  z 
after  superposition  of  membrane  and  bending  stresses  can  easily  be  found. 
The  angles  from  a  reference  direction  (Fig.  2)  to  the  major  membrane  force 
,  to  the  major  plane  of  bending  (M^)  and  to  the  major  total  principal 
stress  a,  at  a  distance  z  are  y,.»  Yk  ®  respectively;  the  angle 

from  to  the  major  bending  plane  (M^)  is  o  and  from  to  the 

major  principal  bending  stress  at  z  is  , 


When  Ah  =  Ah  the  neutral  surface  coincides  with  the  middle  surface; 
h  =  h  -  h  .  'Vhen  Ah  ^  Ah  but  v  =  v*  the  position  of  the  neutral 
surface  is  determined  by: 


h  s 
o 


2  2 

4h^E*/E  +  4hAh  -  (Ah  )  +  (Ah  ) 

_ _ _ u _ o  u 

2(2hE*/E  +  Ah  +  Ah  ) 
o  u 
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«  = 

♦z  =  ■  ®z  ^  ^  °  1 

«|)  =  Yv  -  0  +  90  if  z  <  0 

Z  D  Z  / 


(5) 

(6) 


and  the  positive  quantity  {  ,  identical  for  both  sides,  is  defined  as 


6  =  MAH/NI  s  ^  ° 

The  stress  difference  i  0  at  a  distance  z  is  easily  written  in  terms 

of  and  o,  .  ,  or  with  the  use  of  (2)  and  (7)  in  terms  of  o.  ,  S  and 

t  Jt>l  t 

z 

2 

0^  =  0^  +  2a^aj^^(z/a^)cos2o  +  (Cj^^z/aj,)  (8) 

2 

=  af[l  +  (z«/H)  +  2Cz6/H)cos2a]  (8') 

Z  T 


with  z  <  0  and  i  =  0  for  the  coating  above;  and  z  >  0  ,  i  =  u  ,  below. 

The  angle  <ti  can  be  easily  calculated,  e.g.  from  Mohr’s  circle  in  Figure  4. 

z 

cos2(j)  =  (cos2a  +  6z/H)a  /a  ,  sin2<{'  =  sin2a*a  /o  (9) 

Z  X  Z  Z  X  z 


At  the  free  surfaces  z  =W.  ,  (-H  above,  +H  below)  and  <|>  =<(».«, 

1  ’  o  ’  u  z  i2  ’ 

**’u2^* 


8in24>j^2 


sin2o. 


x/\/iT 


2(«H^/H)cos2a 


> 


(9’) 


cos2<}'jj^2  *  (cos2a  qp.  {H 


2(6H./H)cos2a 


where  the  upper  sign  should  be  taken  with  i  =  o  (coating  above)  and  the 
lower  sign  with  i  =  u  (coating  below).  In  addition,  from  (6) 
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®u2  '  ®o2  ^  *02  "  ♦u2  ^  ^ 

The  proper  quadrant  for  found  from  the  signs  of  their  sine  and 

cosine  in  the  usual  way.  For  2(|»^2  the  regular  angle  must  be  changed  by 
*180°  (Fig.  5). 

It  should  be  noted  that  for  z  >  0  the  positive  bending  stress  differ¬ 
ence  makes  with  the  same  angle  a  as  the  major  bending  plane 

makes  with  N  ,  but  for  z  <  0  the  angle  is  a  +  90  according  to  the  defi¬ 
nition  of  the  stress-difference  as  positive  and  parallel  to  the  major  prin¬ 
cipal  stress.  However  the  same  angle  2a  appears  in  all  expressions  for 
both  coatings,  but  the  proper  sign  is  determined  by  the  value  of  z  ,  ex¬ 
cept  in  expression  (8*)  where  the  coating  dimensions  are  taken  with 

the  proper  sign,  positis«  for  the  coating  below  (i  =  u,  H^>0)  and 
negative  above  (i=o,H^<0). 

Alternatively  the  direction  of  the  bending  stress-difference  could  be 
identified  with  the  plane  of  major  bending  and  its  magnitude  would  be  posi¬ 
tive  for  z  >  0  (when  the  major  stress  is  parallel  to  the  plane  of  bend¬ 
ing)  and  negative  for  z  <  0  (when  the  minor  stress  is  parallel  to  the 

plane  of  bending).  The  final  formulas,  however,  must  be  expressed  in  terms 
of  observed  relative  retardation  or  fringe  order,  which  is  proportional  to 
the  stress-difference,  hence  would  also  be  positive  or  negative.  This 
scheme  would  be  acceptable  for  uniaxial  stress,  as  e.g.  for  a  beam  under 
tension  and  bending  as  in  Fig.  1,  as  well  as  in  ordinary  plane-stress 
photoelasticity  whenever  no  doubt  exists  as  to  the  sign  and  direction  of 
the  stress,  (e.g.  at  the  free  boundary  which  is  the  region  of  main  inter¬ 
est).  In  general,  however,  the  tensorial  variation  cannot  be  escaped  and 
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complications  arise.  If  the  stress  difference  causing  one  fringe  is  defined 
as  positive  when  the  major  principal  stress  has  a  fixed  direction ' and  nega¬ 
tive  when  it  is  perpendicular  to  it,  then  the  same  circular  fringe,  e.g.  in 
the  problem  of  a  shrink  fit  (Fig.  would  be  fringe  +1  at  two  diametral 
points  and  fringe  -1  at  the  diametral  points  of  the  perpendicular  diameter. 
All  other  points  would  have  undetermined  sign,  as  the  major  principal  stress 
has  none  of  the  previous  two  directions.  Obviously  the  first  definition  is 
preferable  and  was  adopted.  The  only  inconvenience  is  the  apparently  dis¬ 
continuous  jump  by  90®  of  the  direction  of  the  stress-difference,  whereas 
under  the  alternative  definition,  it  would  have  only  changed  sign.  In  most 
formulas,  however,  the  same  change  of  sign  could  be  attributed  either  to  an 

always  positive  a  or  n  and  a  changing  cos2a  to  cos2(a  +90)  or  to 

z  z 

a  fixed  cos2a  and  a  sign  change  in  o  . 

z 

Birefringence  in  Inhomogeneous  Stress  Fields 

The  first  study  of  birefringence  in  inhomogeneous  stress  fields  appears 
to  be  the  work  of  Neumann  [10]  who  gave  the  differential  equations  governing 
the  intensity  and  phase  retardation  of  the  components  polarized  along  the 
principal  stresses  for  any  variation  of  azimuth  and  intensity  of  stress- 
difference  along  the  light  path.  Poincare  [11]  in  a  general  study  of  light 
gave  the  most  elegant  representation  of  birefringence  on  a  unit  sphere.  More 
i*ecently  Drucker  and  Mindlin  [12]  solved  the  problem  of  wave  propagation  in  a 
stressed  medium  of  constant  rate  of  rotation  of  the  principal  directions 
per  unit  distance  along  the  ray  (i.e.  as  a  screw  of  constant  pitch)  and  con¬ 
stant  stress-difference ,  which  in  the  absence  of  rotation  would  have  caused 
a  rate  of  change  28  of  phase  difference  per  vinit  distance.  After  neglecting 
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some  terms,  which  are  extremely  small  for  all  practical  photoelastic  mate¬ 
rials  and  rates  of  rotation,  they  reached  explicit  expressions  for  the  am¬ 
plitude  and  the  phase  difference  of  the  components  of  vibration  along  the 
local  principal  stress  axes  at  each  point  in  terms  of  the  dimensionless 
ratios  R  =  fl/8  and  S  =  +(1  +  R^)  .  Each  component  of  vibration  par¬ 

allel  to  a  principal  stress  at  incidence  gives  rise  at  a  distance  x  to 
two  components,  one  parallel  to  the  local  (rotated)  principal  stress  and 
another  transverse  and  (R/S)sin(S0x)  time  smaller.  The  transverse  com¬ 
ponent  vanishes  only  at  retardations  of  integer  wavelengths  (23SX  =  2kiT) 
Accordingly  the  incident  polarization  re-appears  at  integer  wavelengths  at 
the  same  angle  to  the  local  stz*ess  axes  as  at  incidence  and  extinction  may 
be  achieved  with  the  analyzer  likewise  rotated,  but  the  fringes  will  corre 
spond  to  a  birefringence  S  times  larger  than  if  the  stresses  did  not  ro¬ 
tate.  This  may  be  shown  [27]  to  hold  also  for  R  variable  and  for  reflec 
tion  at  integer  plus  one-half  wavelength  retardations  and  return  through 
the  same  path,  as  in  birefringent  coatings.  At  intermediate  retardations 
the  second  component,  if  at  all  significant,  will  produce  an  ellipse  of 
polarization  bearing  little  resemblance  with  the  ordinary  ellipse  in  the 
absence  of  rotation  and  will  render  the  isoclinics  indistinct  and  errone¬ 
ous.  The  second  passage  of  the  light  back  through  the  coating  does  not 
cancel  the  rotational  effects  and  increases  the  complication. 

Drucker's  and  Mindlin's  results  are  a  special  solution  of  Neumann's 
equations  for  a  constant  R  ,  and  though  correct,  show  an  unnecessarily 
large  rotational  effect  because  they  give  the  new  components  of  polariza¬ 
tion  along  the  rotated  stress  axes.  A  part  of  the  correction  is  just  a 
transcription  of  the  ellipse  of  polarization  in  the  new  axes  and  the 
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remalnder  is  the  pure  rotational  effect.  Drucker's  and  Mindlin's  correction 
increases  continuously  with  R  and  S  whereas  the  error  in  the  polariza¬ 
tion  ellipse  in  fixed  coordinates  can  be  shown  to  be  highest  for  R  =  1  and 
to  diminish  to  zero  for  R  tending  to  zero  or  to  infinity  [27].  This  may 
also  explain  why  rotational  effects  were  not  noticeable  in  a  thin  twisted 
tape  subjected  to  light  axial  pull  and  having  a  very  high  rotation  at  a  small 
retardation  [29],  but  became  quite  pronounced  under  a  stronger  axial  pull 
causing  a  reduction  in  R  [27]. 

The  first  application  of  the  effect  of  stress  rotation  was  made  by 
Drucker  [15]  in  a  photoelastic  study  of  plates  subjected  to  bending  at  an 
angle  to  an  initial  frozen-in  tension,  permitting  the  detection  of  the 
otherwise  self-anulling  birefringence  due  to  bending.  The  direction  and 
magnitude  was  known  and  only  the  direction  and  difference  in  principal  bend¬ 
ing  moments,  hence  two  unknowns,  were  sought  from  a  single  photoelastic  ob¬ 
servation.  That  problem  was  simpler  than  the  one  of  the  present  paper  in 
which  4  unknowns  (2  directions,  2  stress-differences)  are  sought  from  2 
photoelastic  observations . 

A  general  study  of  birefringence  in  inhomogeneous  stress  fields  in 
media  with  linear  constitutive  relations  was  done  by  Mindlin  and  Goodman 
[16].  They  have  shown  that  even  in  the  absence  of  rotation  the  ordinary 
photoelastic  law  is  an  approximation  if  the  stress  field  contains  gradi¬ 
ents  along  the  ray,  as  also  shown  by  Mindlin  [14],  or  across  it.  With 
additional  approximations,  similar  to  those  of  Drucker  and  Mindlin  [12], 
they  reached  general  differential  equations  which  for  harmonic  waves  re¬ 
duced  to  Neumann's  equations.  No  explicit  solution  has  been  obtained. 


except  for  R  constant,  but  it  was  shown  that  the  solution  would  mainly 
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depend  on  the  variation  of  R  along  the  light  path.  The  most  systematic 
study  has  been  made  by  Aben  [19,20]  with  matrix  operators. 

It  is  not  always  realized  that  birefringence  with  optical  activity 
(rotating  power)  in  crystals  is  identical  with  birefringence  in  rotating 
stress  fields.  In  the  stress-fields  the  electric  vector  tends  to  remain 
parallel  and  the  principal  directions  rotate,  whereas  in  the  active  crystal 
the  electric  vector  is  rotated  and  the  principal  birefringence  directions 
remain  unchanged.  Accordingly  many  results  obtained  in  physical  optics 
(such  as  the  rotation  of  certain  orthogonal  elliptical  vibrations  without 
change  of  ellipticity)  may  be  directly  applied  to  uniformly  rotating  stress 
fields.  Likewise  the  methods  of  the  Poincare  sphere  [11,17]  and  of  matrix 
operators  [13,17-20]  offer  excellent  means  for  the  visualization  and  calcu¬ 
lation  of  stress  birefringence  effects.  Some  of  these  results  have  fre¬ 
quently  been  re-derived  as  variants  for  special  applications.  An  interesting 
method  consisting  essentially  of  a  plane  projection  of  the  Poincar^  sphere 
was  derived  by  Menges  [21]  and  Kuske  [22].  Wood  [30]  studied  the  vibration 
modes  of  crystal  plates  with  rotating  power.  Mark  [31]  has  calculated  the 

combination  of  bending  and  tension  of  plates  giving  clear  isoclinics.  Stress 
determination  in  rotating  stress  fields  by  scattered  light  was  suggested 
by  Menges  [21]  and  separately  with  a  special  technique  by  Robert  and 
Guillemet  [34,35]. 

The  stress  distribution  in  the  present  problem  is  only  slightly  more 
general  than  in  the  problem  of  the  twisted  tape  [29]:  it  is  a  superpo¬ 
sition  of  a  uniform  tension  on  a  linearly  variable  bending  stress  at  an 
angle  a  which  may  vary  from  0  to  180®,  instead  of  being  always  at  45®. 

On  the  other  hand  the  bending  stress  variation  is  much  smaller  across  the 
coating  than  over  the  whole  coated  shell,  hence  the  rotation  is  frequently 
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small  a  negligibly  so  with  very  thin  coatings.  Certainly  no  rotational  ef¬ 
fect  should  be  expected  at  points  of  the  plate  where  the  principal  direc¬ 
tions  of  bending  and  tension  form  angles  of  0  or  90®  because  the  linear 
distribution  of  stresses  along  the  light  path  do  not  give  rise  to  measur¬ 
able  discrepancies  firom  the  simple  photoelastic  law  [14,  16],  Further¬ 
more,  when  either  stress  is  much  stronger  than  the  other,  the  rotation  is 
very  small  and  without  effect  on  the  relative  retardation.  Some  indica¬ 
tion  of  the  probable  rotational  effect  on  the  birefringence  may  be  obtained 
from  the  average  value  of  R  across  the  coating  thickness  (calculated  as 
the  ratio  of  the  rotation  to  half  the  retardation  for  single  passage  through 
the  coating  with  uniform  stress  equal  to  the  stress  at  mid-thickness). 

Figures  3a,  3b  show  graphs  of  vs.  the  angle  a  for  membrane  and  bend¬ 

ing  fringe  orders  n^  and  nj^  respectively  in  0.108  in.  (Fig.  3a)  or  0.057in, 
thick  coatings  (Fig.  3b)  on  a  0.250  in.  thick  shell.  Obviously  significant 
rotational  effects  should  be  expected  when  n^  and  n^^  are  small  and  about 
equal  and  simultaneously  the  angle  o  is  between  about  70®  and  110®  but  not 
very  close  to  90®,  as  then  no  appreciable  error  is  expected  (all  the  rotation 
then  occurs  in  the  region  where  the  stress  is  negligible ,  hence  R  is  very 
large  and  the  effect  small). 

The  rotation  and  change  of  magnitude  of  the  total  principal  stress 
throTigh  the  coating  thickness  can  be  easily  visualized  on  a  Mohr  diagram. 

As  the  photoelastic  effect  depends  on  principal  stress  differences,  it  is 

permissible  to  consider  both  membrane  and  bending  stresses  as  equivalent 

*t  t  1  b  Id  X 

states  of  pure  shear,  o,  =  -a„  =  w  o.  and  o,  =  -a-  =  tt  Oi,  • 

1  2  2  t  Iz  2z  2  DZ 

anticlastic  bending  this  is  exact;  for  other  problems  the  true  state  will 
differ  by  an  isotropic  stress,  which  anyway  cannot  be  determined  from 
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simple  photoelastic  observations  at  normal  incidence  at  a  point.  These 
assumptions  may  shift  the  origin  of  Mohr's  circle  but  do  not  affect  its 
diameter  and  they  greatly  simplify  the  superposition  of  the  stress-differ¬ 
ences.  Figure  4  (top)  shows  the  separate  membrane  and  bending  stress 
states  and  their  superposition  at  a  general  point  in  the  coating  below 
(z  >0)  and  in  a  symmetric  point  (  z' <  0)  in  the  coating  above  as  seen 
from  some  point  z  »  0  .  The  membrane  stress  is  identical  at  the 

two  points,  but  the  bending  stress-difference  ,  positive  by  defini¬ 
tion,  is  at  an  angle  a  to  for  z  >  0  and  at  a  +  90  for  ■z*<  0  . 

The  superposition  is  carried  out  in  Mohr's  diagram  by  drawing  first  the 
circle  for  >  0  (Fig.  4),  determining  the  points  A  and  B  on  the 

circumference  at  angles  2a  and  2(o  +  90®)  and  adding  the  suitable 
bending  stresses  i  and  -  ^  to  find  the  points  A^,  (coat¬ 

ing  below)  or  A^,  B^  (coating  above),  which  represent  the  total  second¬ 
ary  stresses  in  the  principal  bending  directions.  Mohr's  circle  drawn  on 
A  B  (or  A  B  )  as  diameter  represents  the  total  state  of  stress  at  z 

(or  z*)  and  these  diameters  represent  the  corresponding  total  stress- 

difference.  The  angle  ^  from  the  major  principal  total  stress  a,  to 

z  xz 

the  major  bending  stress  is  also  shown. 

The  variations  of  stress  through  the  coating  thickness  can  be  found 
on  a  similar  diagram  simply  by  varying  the  magnitude  of  the  bending  stress 
^  ,  i.e.  the  lengths  of  AA^  and  BB^  ,  etc.  For  example.  Figure  5 

shows  the  superposition  of  membrane  and  bending  stress  differences  at  the 
free  surface  (diameter  ^y2®u2^  metal-coating  interface  (diame¬ 
ter  coating  below,  and  similarly  for  the  coating  above 

(diameters  A  _B  .  and  A  ,B  , 
o2  o2  ol  ol 


respectively).  The  corresponding  Mohr's 
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circles  need  not  be  drawn.  The  angles  major  total 

stress  at  the  free  surface  and  at  the  interface  of  the  coating  below  with 
the  major  bending  stress  are  easily  found,  and  also  the  corresponding 
angles  and  4)^^^  in  the  coating  above.  The  total  rotation  ~ 

across  the  thickness  of  the  coating  above  is  much  larger  than  the  rotation 
*^o2  ~  ^ol  orientation  a  of  in  relation  with  ,  i.e. 

the  position  of  point  A  ,  greatly  influences  the  total  rotation  through 
the  coating,  evenfbr  constant  and  . 

The  points  P^,  represent  the  stress  state  at  a  distance  z  >  0  . 
Their  position  along  the  segments  ®xil^u2  proportion¬ 

ally  with  the  bending  stress-difference ,  hence  proportionally  with  z  . 
Twice  the  length  OP^  gives  the  corresponding  total  principal  stress- 
diffei?ence  and  half  the  angle  DOP^  gives  the  angle  of  to  the 

constant  direction  of  .  Obviously  point  gives  the  lowest  stress-^ 

difference  and  the  highest  rate  of  rotation  per  unit  change  of  z  ,  hence 
gives  the  highest  ratio  R  .  When  this  point  lies  in  the  coating  and  the 
distance  from  0  to  is  small,  the  total  rotation  through  the  coat¬ 

ing  should  be  large  and  the  stress  small,  hence  the  effects  on  polariza¬ 
tion  could  be  significant.  In  the  coating  below  this  situation  arises  when 
a  approaches  90®  (point  A  close  to  C  )  and  simultaneously  the  bending 
stress  at  the  coating  mid-thickness  is  about  equal  to  .  For  the  coat¬ 
ing  above  a  should  approach  zero.  This  was  also  the  conclusion  reached 
from  the  graphs  of  o  (Fig.  3a,  b),  but  it  must  be  remembered 

that  angles  equal  to  or  very  close  to  90°  or  0°  will  cause  no  rotational 


errors . 
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It  is  clear  that  the  region  of  highest  rotation  may  lie  outside  both 
coatings  and  then  the  effects  of  rotation  on  fringe  order  and  on  isoclinics 
may  be  very  small.  Even  when  this  region  lies  within  a  coating,  the  effect 
may  not  be  severe,  as  e.g.  in  the  early  tests  with  a  twisted  tape.  These 
conclusions  spurred  the  hope  that  a  practical  and  sufficiently  accurate 
method  for  finding  the  stresses  in  shells  might  be  developed  even  though 
the  exact  effects  of  rotation  are  not  taken  under  consideration.  The  pro¬ 
posed  method  is  important  and  useful  not  so  much  for  finding  the  birefrin¬ 
gence  once  the  stresses  are  known,  which  can  be  easily  foimd  by  integrating 

✓ 

Neumann’s  equations  by  steps  (or  using  Poincare's  sphere  or  matrix  opera¬ 
tors  ) ,  but  for  the  inverse  problem  of  finding  the  stresses  from  the  bire¬ 
fringence  on  the  two  coatings.  The  integrated  effects  of  birefringence  and 
rotation  in  the  two  coatings  cannot  be  solved  for  the  stresses ,  whereas  the 
proposed  simplified  algebraic  expressions  can.  The  accuracy  of  the  proposed 
method  has  generally  been  found  good,  though  somewhat  less  so  with  a  few 
combinations  of  M  ,  N,  a  causing  a  high  rotation  in  one  coating. 


Calculation  of  the  Birefringence 


The  fringe  order,  if  independent  of  rotation,  would  be  found  as  the  in¬ 
tegrated  photoelastic  effect  during  the  forward  and  backward  passage  through 
each  coating  of  thickness  h^ 


n.  = 
1 


2 

C 

a 


fH. 

1 

o  dz 
z 


(11) 


h. 

1 


-16- 


with  i  =  o  for  the  coating  above  and  i  =  u  below.  Substitution  of  a 

z 

as  in  (8*),  integration  and  substitution  of  fringe  orders  for  stresses 
from  (3),  gives 


n.  =  n  . (H/26Ah,)|A.  -  B.  +  sin^2a]l  n(A./B.)|  (lla,b) 

iti  I'll  ii' 

and  in  terms  of  an  always  positive  H.  (i.e.  H^>0,  H^>0) 

A^  =  ( 6H^/H  +  cos2a)  ^  +  sin^2a 

B.  =  (6h./H  T  cos2o)  \/(6h./H  +  cos2a) t  sin^2a 

Equation  (11a)  for  i  =  o  and  the  upper  of  the  double  signs  gives  the  in¬ 
tegrated  fringe  order  n^  for  the  coating  above;  Eq.  (11b)  with  i  =  u  and 

the  lower  signs  gives  the  fringe  order  below.  The  factor  (n^.H/6Ah^)  is 
the  same  for  both  coatings  and  equal  to  NH/A6C^  which  may  be  substituted 
if  n^  must  be  expressed  in  terms  of  N  instead  of  n^  .  For  identical 

coatings  h./H  =  h/H  ,  H./H  =  1  and  n. .  =  n^  ,  but  n  ^  n  .  These 
1  *1  tit  o  u 

equations  are  more  complicated  than  the  one  given  by  Drucker  [15],  as  they 
contain  the  three  basic  parameters  n^  (or  N),  6  ,  and  a  as  well  as 
the  ratio  h/H  . 

Inverse  Solution  for  Membrane  and  Bending  Stress-Differences 

In  the  inverse  solution  it  is  required  to  determine  n^^,  n^,  o  and 

,  or  the  equivalent  6  ,  n^,  a,  Yj^  ,  from  the  experimental  measures  of 

n  ,  n  ,  0  and  6  _  .  However,  once  6  and  a  are  foxmd  y,  is  easily 

o  u  o2  u2  '  b 

determined  from  (9')  and  (6).  Essentially  the  problem  then  consists  of 

finding  the  three  quantities  6,  n.  and  a  in  terms  of  n  ,  n  and 

t  o’  u 
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6^2  “  ®q2  requires  three  independent  equations  in  these  variables.  If 

equations  (lla,b)  are  used  the  solution  will  have  to  be  numerical.  A  very 
simple  explicit  solution  can  be  obtained  with  a  further  approximation  which 
is  justified  by  the  test  results.  For  thin  coatings  the  principal  stress 
difference  is  assimed  linear  across  the  coating  thickness.  With  this  assump¬ 
tion  the  observed  fringe  order  n^  may  be  taken  equal  to  the  fringe  order  of 
a  coating  with  a  uniform  stress  equal  to  the  real  stress  at  mid-thickness 
z  =  a^  .  Two  equations  are  obtained  from  expression  (7)  for  z  =  a^  after 
multiplication  throughout  by  2Ah/C^  to  transform  stresses  into  fringe 
orders : 


n?^ 

1 


2  2  — 

n^.  =  n^.  + 
mi  ti  ^ 


2n^.n,  .  cos2a  + 
ti  bi 


n 


bi 


(12) 


with  i  =  o  and  the  upper  sign  for  the  coating  above ,  or  i  =  u  and  the 
lower  sign  for  the  coating  below. 

For  identical  coatings,  n^^  =  n^  ,  nj^^  =  n^^  and  a^  =  a  ,  but 
n^  /  n^  .  Addition  and  stibtraction  and  substitution  of  (7)  gives 

n^  +  n^  =  2n^(l  +  6^a^/H^) 

u  o  t 

► 

n^  -  n^  =  4n£( 6a/H)cos2a 

u  o  t 


(13) 


Division  of  the  first  by  the  second  Eq.  (13)  gives 


where 


co=2a  =  m 

25a/H 

m  =  (n^  -  n^)/(n^  +  n^) 
u  o  o  u 


(14) 


(15) 
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The  third  independent  equation  is  obtained  from  the  calculation  of 
cos2(6^2  “  ®o2^  after  the  substitution  (10),  expansion  and  substitution 
of  (9')  for  =  H  ,  or  directly  from  a  Mohr  diagram  as  in  Fig.  4. 


cos(e^^2  “  ®o2^  =-cos((|»^2  ”  =  (1  -  5^)/  \/(l  +  6^)^  -  452cos22a  (16) 


Substitution  of  cos2a  from  (14)  and  solution  for  6^  gives: 


Cl  -  l-m2  cos2(0  -  -  e  )]2_  A 
^2  -  _ o2 _ 

1  -  (l-m2a2/H2)cos 22(6^2  "  ®o2^ 


(17) 


where  A  is  the  quantity 


A  =  2cos2(0  _  -  0  - 
u2  o2 


l-m^-m^sin22( 


®u2  “ 


) 


(18) 


The  value  of  A  is  zero  for  the  following  cases: 


I.  m=0,  i.e.  n  =n  ,  which  occurs  when 
’  o  u  ’ 

a)  n,  =  0  and  therefore  n  =  n  =  n.  and  0  „  -  0  „  =  0 
b  o  u  t  u2  o2 


b)  n.  =  0  and  therefore  n  =n  =n.  2(0--0-)=  ±180® 

t  o  u  b  u2  o2 

c)  n.  =  n,H/a  .and  it  can  then  be  shown  that  2(0  _  -  0  -)  =  *90® 

ZD'  U2  02 

for  any  value  of  a  . 


II.  m  =  l,  i.e.  n^  =  0,  which  occurs  when  n^  =  nj^  and  a  =  0  . 

Then  the  angles  0  „  and  therefore  2(0  „  -  0  -)  cannot  be 
o2  u2  o2 

experimentally  determined,  but  for  a  =  0  Eq.  (16)  shows  that 
sln2(e^j  -  =  0  . 
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III.  m  =  -1  ,  i.e.  =  0  ,  which  occurs  when  s  and  a  -  90®  . 
Again  Eq.  (16)  shows  that  sin2(0^2  ~  ~  ^  ‘ 

The  value  of  A  is  also  zero  for  intermediate  values  of  ,  -1  <  m  <  0 

and  0  <  m  <  1  ,  when  the  coating  is  infinitesimal.  For  thin  coatings 
(Ah  «  H)  A  is  small  and  has  been  neglected  in  all  present  calculations. 

The  value  of  A  may  serve  as  a  check  of  the  experimental  data.  With 
the  correct  values  of  m  and  ~  ®q2  quantity  under  the  second 

square  root  in  A  should  be  positive,  since  a  real  6  always  exists  (in¬ 
cluding  6  -*■  <»  for  n^  =  0). 

1  -  m^  -  m2sin22(0  -  0  -)Ah2/4aH  i  0  (19) 

Experimental  data  not  fulfilling  (19)  would  indicate  measurement  errors  or 
appreciable  rotational  effects.  One  could  then  accept  a  further  approxima¬ 
tion,  either  by  "conditioning"  the  data  to  fulfill  (19),  or  by  neglecting 
the  quantity  A  .  In  the  present  tests  A  was  neglected,  but  only  after 
checking  for  the  fulfillment  of  (19)  by  each  set  of  measiirements .  Violation 
of  (19)  could  possibly  occur  at  angles  a  close  90®  for  the  coating  below, 
or  0®  above  and  for  ^  '^t  ^  Rj^H/a  ,  when  m  approaches  1  but  is  not  so 
close  as  to  make  A  almost  zero  (as  in  II  and  III  above).  These  are  also 
the  conditions  giving  the  strongest  rotation  so  that  errors  would  not  be  un¬ 
expected. 

After  6  is  found  from  (17),  the  value  of  n^  is  found  from  the  first 

of  (13)  and  of  n,  from  (7) 
b 

nl  =  (n2  +  n2)/(l  +  62a2/H2)  (13) 

T  z  U  O 


n,  =  n. 6H/a 
D  t 


(7) 
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Finally  the  value  of  a  is  found  from  (14).  Calling  2a  the  1st  quadrant 
solution  of  (14)  when  cos2o>0  or  m  >  0  or  the  2nd  quadrant  solution 
when  cos2a  <0  or  m  <  0  ,  the  solution  for  2a  is 


2a  =  2a  when 

“  ‘  ^<®u2  -  ®o2'  ‘  ’ 

or 

2a  =  -2a  when 

<  2(6^2  -  «o2'  *  ” 

or 

0  <  2(4.„2 

-t.2> 

The  computed  values  of  n^  ,  and  a  give  the  magnitude  of  the  membrane 
and  bending  "stress-difference”  and  their  relative  angle.  Their  absolute 
directions  Yjj  can  easily  be  found  from  (9'),  (6)  and  (5),  as  already 

mentioned.  Instead  of  either  of  (9’)  it  is  easier  to  use  their  quotient  to 
find  ♦^2* 


tan2+^2  “  sin2a/(cos2a  +  6) 


(20) 


with  the  usual  quadrant  selection  according  to  the  signs  of  numerator  and  de¬ 
nominator. 


Test  Methods 

The  tests  had  two  purposes:  to  check  the  accuracy  of  the  direct  formulas 
(8')  and  (11)  or  (12)  giving  the  total  birefringence  in  terms  of  the  applied 
membrane  and  bending  stresses  and,  if  these  proved  accurate  enough,  to  check 
the  inverse  formulas  (15)-(17)  and  (10')  giving  the  applied  membrane  and  bend¬ 
ing  stresses  (or  corresponding  birefringence  and  their  relative  direction)  in 
terms  of  the  observed  fringe  orders  and  principal  directions  on  the  two  coat¬ 
ings.  Two  types  of  tests  were  made,  both  with  doubly  coated  square  plates  sub¬ 
jected  to  anticlastic  bending  and  to  in-plane  tension  of  variable  direction. 
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In  both  series  of  tests  anticlastlc  bending  was  produced  by  transverse 
forces  Q  applied  at  a  pair  of  diagonally  opposed  comers  and  -Q  at  the 
other  pair,  which  caused  principal  bending  moments  =  -M2  "  ^  Q  or 
M  =  Q  ,  with  principal  directions  parallel  to  the  plate  diagonals. 

First  Series  of  Tests.  In  the  first  series  in-plane  tension  was  produced 
by  forces  P  evenly  distributed  through  a  system  of  beams  and  levers  at 
16  points  of  the  half-perimeter  of  the  plate,  8  on  each  side,  at  an  angle  a 
to  the  diagonals  which  could  be  varied  between  45°  and  -45°  (Figs.  6,  7). 

Anticlastic  bending  deformed  the  plate  into  a  h3rperbolic  paraboloid. 

To  cause  only  membrane  stress  the  16  loads  on  the  perimeter  should  be  tan¬ 
gential  to  the  deformed  plate.  This  could  be  approximated  by  an  out-of- 
plana  tilt  of  the  system  of  levers  about  the  midpoints  H  (Fig.  6),  ad¬ 
justed  after  several  trials  to  the  point  where  a  change  of  P  from  zero 
to  its  highest  value  would  cause  no  variation  of  the  already  applied  Q  . 

The  adjustment  was  generally  difficult,  lengthy  and  imperfect,  except  when 
a  was  equal  to  0°  or  45°.  As  a  consequence  an  tinknown  small  stress  system 
was  introduced  and  caused  errors.  Other  sources  of  error  were  the  imperfect 
centering  of  the  pins  about  the  mid-plane  of  the  plate  and  of  the  levers , 
and  the  unavoidable  location  of  the  pins  in  holes  lying  a  short  distance 
inside  the  edges  of  the  plate.  Halfway  between  the  holes  the  edge  was 
slotted  to  the  depth  of  the  centerline  of  the  holes,  which  thus  determined 
the  size  of  the  square  subjected  to  tension.  On  the  other  hand,  to  achieve 
an  exact  anticlastic  loading  over  the  whole  plate ,  the  loads  Q  had  to  be 
applied  precisely  at  the  comers  of  the  original  square ,  which  was  larger 
than  the  one  subjected  to  tension.  A  compromise  was  reached  by  slotting 
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only  the  plate  and  not  the  birefringent  coatings  so  that  at  the  slots  the 
plate  was  weak  in  tension  but  could  take  a  good  part  of  the  bending.  There¬ 
fore,  to  a  first  approximation,  the  tension  was  applied  to  an  8"  square  and 
the  bending  to  a  9”  square. 

The  plate  was  0.248  in.  thick  24ST  aluminum,  and  the  coatings  were 

1” 

selected  relatively  stress-free  ^  thick  sheets  of  CR39.  The  small  sen¬ 
sitivity  and  considerable  thickness  permitted  an  assessment  of  the  present 
method  under  unfavorable  conditions. 

The  values  of  the  various  quantities  were; 


2h 

=  0.250  in 

H 

o 

0.3705  in 

E  = 

2.86  X  105  psi 

Ah 

=  0.245  in 

H 

0.3955  in 

V  = 

0.40 

o 

u 

Ah 

=  0.271  in 

a 

= 

0.248  in 

A  = 

9.90  in^/in 

u 

o 

h 

=  0.1255  in 

a 

= 

0.260  in 

I  = 

0.853  in*'/in 

o 

u 

+ 

h 

u 

=  0.1245  in 

E* 

1.02  X  10^  psi 

C  = 
a 

78  lb/in(in.  fringe) 

.2H 

=  0.766  in 

V* 

= 

0.33 

The  inequality  of  the  two  coatings  had  a  very  small  influence  on  the 
position  of  the  centroid  (0.001"  off  mid-thickness)  and  an  altogether  neg¬ 
ligible  effect  on  the  effective  moment  of  inertia.  All  calculations  were 
made  with  the  expressions  for  assjrmetric  coatings. 

These  first  tests  served  to  verify  expressions  (11)  or  (13')  for  total 
birefringence  and  (8')  for  its  principal  direction.  Observations  were  made 
only  on  the  coating  below.  Plane  polarized  light  was  used  at  the  azimuth 
giving  the  brightest  fringe  pattern.  Rather  than  apply  arbitrary  values  of 


Determined  by  a  test  in  pure  anticlastic  bending. 
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N  and  M  (or  n.  and  n,  )  and  2a  and  then  check  the  observed  vs. 

tU  DU 

the  calculated  values  of  n^  and  6^2  >  ^  was  applied  first  and  then  N 

was  increased,  at  a  constant  direction,  till  the  first  fringe  order  was 
observed.  For  i  =  u  ,  the  lower  sign  (+)  and  n^  =  1  ,  Equation  (12)  rep¬ 
resents  a  family  of  central  ellipses  with  semi-axes  along  the  bisectors  of 
the  coordinate  axes  nj^,  n^  equal  to  0.707seca  and  0.707coseca  ,  and 
passing  through  the  points  (n^^  =  ±1  ,  n^  =  0)  and  (nj^  =  0  ,  n^  =  *1). 
Figure  8  shows  these  ellipses  and  the  corresponding  experimental  points 
(i.e.  the  values  n^,  which  in  the  test  caused  one  fringe  at  angles  a 
equal  to  0®,  15®,  30®,  45®  and  90®). 

As  the  principal  directions  could  not  be  easily  detected  when  the 
whole  plate  was  covered  with  the  dark  first  fringe,  the  loads  were  re¬ 
duced  by  ^  and  the  crossed  polaroids  were  rotated  to  give  minimum  light 
intensity.  In  irrotational  photoelasticity  this  would  have  given  the  true 
isoclinic  at  y  as  well  as  1  fringe. 

At  angles  a  equal  to  60®  and  75®  the  tension  needed  to  produce  one 
fringe  exceeded  the  capacity  of  the  machine  even  when  the  bending  was 
light.  Tests  were  then  made  with  polarizer  and  analyzer  parallel  (light 
field)  and  membrane  and  bending  intensities  adjusted  to  give  a  black  fringe, 
which  in  ordinary  irrotational  photoelasticity  would  correspond  to  a  fringe 
order  of  ^  .  Of  course  the  rotational  effect  at  a  retardation  of  —  wave¬ 
length  is  stronger  than  at  one  wavelength,  and  must  have  influenced  the 
fringe  order  as  well  as  the  direction.  In  addition  the  90®  rotation  of  the 
analyzer  introduced  further  complications,  but  it  was  still  desirable  to 
test  the  suitability  of  this  approximate  method  under  such  conditions. 
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Agreement  between  analytical  and  experimental  fringe  orders  was  excel¬ 
lent  at  0  =  0°  and  reasonably  good  at  o  =  15°  ,  30°  and  in  part  at  45° , 
with  loads  producing  a  retardation  of  one  fringe  (n  =  1)»  but  not  so  good 
at  o  =  60°  and  75°  and  with  the  lower  loads  producing  only  half  a  fringe 
(n  =  ^).  A  significant  part  of  the  errors  was  caused  by  pin-loading  in 
tension,  as  already  discussed  and  indicated  by  the  xmcertainty  in  the  load 
needed  to  cause  one  fringe  even  in  simple  tension  (Fig.  8,  points  on  verti¬ 
cal  axis ) ,  when  the  plate  had  no  curvature .  The  errors  were  increased  by 
the  bending  curvature,  by  appreciable  initial  stresses  and  by  viscoelastic 
effects,  especially  after  prior  loadings.  The  principal  directions,  how¬ 
ever,  showed  such  a  high  discrepancy  that  no  attempt  was  made  to  solve  in¬ 
versely  for  the  applied  membrane  and  bending  stresses  in  terms  of  the  ob¬ 
served  birefringence  on  the  two  coatings. 

Second  Series  of  Tests 

a.  The  Plate.  The  purpose  of  the  second  series  of  tests  was  to  increase 
the  accuracy  and  check  the  inverse  solution  for  membrane  and  bending  loads 
in  terms  of  observed  birefringence  on  the  two  faces.  The  tests  were  there¬ 
fore  designed  so  as  to  remove  the  uncertainty  of  pin-loading  and  to  reduce 
the  errors  due  to  the  excessive  thickness,  to  the  low  sensitivity  and  to 
the  viscoelastic  behavior  of  the  previously  used  coatings.  To  avoid  pin¬ 
loading  or  clamping  and  the  resulting  transverse  constraint,  it  was  decided 
to  simulate  the  membrane  stress  by  a  uniform  residual  tension  imposed  on 
the  coatings  before  cementing  them  to  the  plate.  Anticlastic  bending  would 
be  applied  as  before.  The  difficulty  of  uniformly  stressing  a  9  x  9  in. 
plate,  however,  would  not  be  overcome  but  only  transferred  from  the  plate  to 
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the  coating.  Furthermore  pairs  of  identical  coatings  would  have  to  be 
stressed  to  various  intensities  and  be  cemented  with  their  principal  di¬ 
rection  at  various  angles  to  the  diagonals.  It  was  then  recognized  that 
the  strived-at  uniformity  of  stress  over  the  whole  plate  was  not  essen¬ 
tial,  since  each  test  required  an  observation  at  only  one  point  on  each 
coating.  Pairs  of  coatings  with  identical  residual  stress  fields,  but 
variable  so  as  to  include  all  required  magnitudes  and  directions,  would 
permit  the  use  of  a  single  carefully  coated  plate  for  all  tests.  Each 
test  would  again  consist  of  observations  at  two  corresponding  points, 
one  on  each  coating,  having  the  desired  magnitude  and  direction  of  mem¬ 
brane  stress.  Although  in  the  present  tests  the  residual  stress  fields 
of  the  two  faces  were  identical  to  a  high  degree  of  accuracy,  this  re¬ 
quirement  could  have  been  relaxed,  since  anticlastic  bending  was  uniform 
over  a  large  central  area  of  the  plate.  With  the  dimensions  and  loads 
used  the  errors  from  edge  effect,  membrane  action  and  transverse  shear 
[323  should  be  less  than  0.002  in.  a  central  5  in.  diameter  circle.  Ob¬ 
servations  could  then  be  made  at  any  two  points ,  one  on  each  face ,  having 
identical  membrane  stress ,  even  if  they  would  not  lie  on  the  same  normal 
to  the  mid-plane. 

Several  methods  were  tried  for  casting  the  araldite  sheets  used  as 
coatings.  Most  methods  caused  residual  stresses,  cracking,  or  an  uneven 
thickness.  The  final  and  most  successful  method,  perfected  after  many 
trials,  used  specially  machined  13  x  15  x  3  in.  aluminum  molding  plates 
ground  flat  within  *0.0005  in.  and  then  polished  and  chrome-plated  to  a 
surface  smoothness  of  10  ^  in.  Two  matching  pairs  of  such  plate  molds 


were  used.  A  thin  coat  of  a  mixture  of  10%  Dow  Silicone  20  mold  release 
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grease,  10%  toluene  and  80%  isopropanol  was  applied  and  wiped  off  after  dry¬ 
ing.  The  mold  plates  were  spaced  by  0.5  in.  wide  annealed  strips  of  cellu- 

_ lose  acetate  millejl  to  a_iinifomu  thickness  5nd^«ej5eL  hel<Ljagainst  each  other _ 

by  10  screws  on  the  perimeter.  The  spacers  were  placed  first  on  three  sides 
leaving  the  top  open  for  pouring  the  hot  araldite ,  after  which  the  fourth 
spacer  was  put  in  place  and  the  screws  were  uniformly  and  suitably  tightened. 
The  function  of  the  plastic  spacers  was  to  yield  in  coiiiprearsion  as  they 
softened  and  permit  a  reduction  of  thickness  matching  more  or  less  the  volume 
contraction  of  araldite  during  setting.  In  contrast  with  metal  spacers  which 
had  caused  surface  separation  during  casting  as  well  as  high  residual  stresses 
and  cracks,  the  plastic  spacers  gave  excellent  sheets  of  uniform  thickness 
having  negligible  or  annealable  stresses.  After  several  castings  it  was  pos¬ 
sible  to  select  pairs  of  9  in.  square  sheets  with  thickness  variations  no 
larger  than  ±0.0005  in.  within  each  sheet  or  from  each  other,  which  was  con¬ 
siderably  better  than  with  any  commercially  available  product.  Sheets  of 
thickness  0.108,  and  0.057  in.  were  chosen. 

Araldite  6060  resin  with  30%  by  weight  of  phthalic  anhydride  was  mij^ed 
at  140°C,  filtered  through  fiberglass  cloth  and  cast  in  hot  molds  at  120®  ± 

5®C  and  cured  at  115®C  for  8  hours.  The  cast  sheets  were  removed  cold  and 
were  laid  on  an  open  paper-covered  mold  and  weighted  with  a  thick  rubber 
sheet  for  annealing  at  140®C  for  3  hours  followed  by  cooling  at  the  rate  of 
50/h  to  room  temperature. 

The  residual  stress  field  was  produced  by  ’’shrink-fitting,”  a  1  in. 
disc  in  a  narrower  hole  bored  centrally  in  each  sheet.  The  disc  was  ma¬ 
chined  from  the  border  of  the  plate  (initially  13  x  15  in;  finally  9  x  9  in.) 
to  which  it  was  fitted.  Both  turning  of  the  disc  and  boring  of  the  hole  had 


-27- 


to  be  done  with  the  greatest  care  on  the  most  accurate  machines  to  achieve 
perfect  roundness  and  to  prevent  chatter-marks  and  uneven  fringe  patterns. 
The  interference  (0.005  to  0.010  in.)  was  such  as  to  produce  about  8  con¬ 
centric  fringes  in  reflected  light  in  the  0.108  in.  coating  and  about  5  in 
the  0.057  in.  coating.  For  the  two  thicker  coatings  the  discs  were  cooled 
in  liquid  nitrogen  and  were  quickly  fitted  in  the  plates  which  had  been 
heated  to  about  45®C.  The  same  procedure  caused  conical  buckling  in  the 
thinnest  coating.  Accordingly  part  of  the  fringe  pattern  was  ^ozen  jin  the 
sheet  by  forcing  in  the  hole  of  the  hot  clamped  sheet  a  tapered  round  bar 
ending  in  a  cylindrical  portion  of  the  desired  diameter.  The  frozen-in 
pattern  was  increased  by  a  residual  stress  pattern  by  shrink-fitting  the 
disc.  This  ensured  also  sufficient  pressure  between  disc  and  sheet  to 
prevent  separation  during  unfavorable  bending. 

The  pairs  of  sheets  were  finally  cemented  on  both  faces  of  the  0.250 
in.  thick  plates  of  24ST  aluminum.  These  plates  were  square  with  9  in. 
sides  and  had  small  protrusions  at  the  comers  in  the  directions  of  the 
diagonals  permitting  the  transverse  load  to  be  applied  exactly  at  the  cor¬ 
ners  of  the  square.  Transverse  pins  were  fitted  in  these  protrusions, 
centered  on  the  comers  and  protruding  by  0.125  on  each  side,  so  as  to 
support  spherical  ball  bearing  through  which  the  transverse  loads  were 
applied  (Fig.  9).  Cementing  was  done  with  cold  setting  araldite  in  a 
thick  layer  on  which  the  well  centered  9  in.  square  sheets  were  placed 
and  then  uniformly  weighted  so  as  to  expell  the  excess  glue.  Curing 
lasted  for  a  week  at  room  temperature.  Lack  of  miformity  of  the  glue 
line  unfortunately  reduced  the  high  accuracy  of  the  sheet  thickness.  The 
coatings  had  a  final  thickness  of  0.108  *  0.0015  and  0.0565  ±  0.0015  in. 
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b.  Measurements.  Observations  were  made  at  normal  incidence  with  the  help 
of  a  set  of  3  full  mirrors  and  a  partially  reflecting  front  surface  mirror 
placed  in  the  field  of  a  regular  transmission  polariscope  as  shown  in  the 
diagram  of  Fig.  10  and  the  photograph  of  Fig.  11.  Plane  polarized  light 
with  a  fixed  plane  of  vibration  perpendicular  to  the  common  plane  of  inci¬ 
dence  of  all  mirrors  was  used  throughout  so  as  to  avoid  the  depolarization 
from  reflexion  at  other  azimuths.  The  various  isoclinics  were  obtained  by 
rotating  the  plate  in  its  plane.  The  angular  position  was  read  on  a  divided 
circle  below  the  plate  (Fig.  9).  Fringe  measurements  were  made  to  the 
nearest  0.1  fringe  at  the  plate  azimuth  making  brightest  the  neighborhood  of 
the  point  under  observation.  ”Isoclinics"  were  determined  to  the  nearest  1®, 
but  sometimes  contained  an  uncertainty  of  *5®  or  more.  The  accuracy  of  the 
measurements  was  purpose dly  kept  at  the  level  of  ordinary  visual  observations 
without  compensator  or  photometer. 

As  already  mentioned  it  was  found  more  convenient  to  express  and  cal¬ 
culate  bending  moment  and  membrane  force  in  terms  of  the  fringe  order  in  the 
coatings,  but  with  expressions  (3)  all  results  may  be  rewritten  in  terms  of 
M  and  N  .  The  fringe  order  from  bending  alone  could  be  easily  detected  in 
the  central  expansion-fitted  disc  which  initially  was  in  an  isotropic  state 
of  compression,  hence  had  a  zero  fringe.  A  more  accurate  method  was  to  meas¬ 
ure  the  radial  displacement  of  the  "membrane  stress"  fringes  along  the  plate 
diagonals,  where  membrane  and  bending  directions  are  parallel.  Displacement 
of  the  integer  membrane  fringe  n^  to  the  initial  position  of  the  fringe 
n^  -  1  along  one  diagonal  and  to  n^  +  1  along  the  other  indicates  a  value 
of  equal  to  1  ;  displacement  of  n^  to  n^  -  2  and  n^  t  2  a  value 

nj^  =  2  ,  etc.  as  shown  in  Figures  12a-c  of  the  lower  coating  (z  >  0) 
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specially  taken  in  circularly  polarized  light.  For  measurements  on  the 
opposite  side  (z  <  0)  the  plate  was  turned  over  and  subjected  to  the  same 
bending,  as  if  it  had  never  been  \inloaded.  To  achieve  identical  bending 
in  both  positions  reliance  was  placed  not  on  the  intensity  of  the  comer 
loads,  measured  by  a  dynamometer,  but  on  the  plate  deflection  indicated 
by  a  special  very  sensitive  deflectometer  spanning  one  diagonal  (Fig.  9). 

In  other  words  bending  was  controlled  by  fixed  boundary  displacements  in¬ 
stead  of  fixed  loads.  Consequently  any  photo-viscoelastic  fringe  shift 
would  be  a  relaxation  effect,  which  has  been  fomd  [29]  to  be  appreciably 
smaller  than  the  corresponding  creep  under  constant  applied  loads. 

At  any  point  of  the  plate  the  major  membrane  stress  is  circinnferen- 
tial  (shrink  fit  causes  a  hoop  tension  and  a  radial  compression  outside 
the  central  disc)  and  the  major  principal  bending  plane  is  everywhere  (ex¬ 
cept  near  the  edges)  parallel  to  the  diagonal  along  which  the  initially 
circular  fringes  contract.  For  convenience  this  major  diagonal  on  the 
lower  coating  was  chosen  as  the  reference  direction,  hence  in  the  present 

tests  Yjj  =  0  ,  o  =  -  and  =  -6^2  +  ’  ♦u2  =  '■®u2 

ure  13  shows  corresponding  points  in  the  lower  and  upper  coatings  having 
a  membrane  "stress"  of  one  fringe  with  principal  direction  at  an  angle 
q  =  30®  to  the  major  plane  of  bending  (diagonal  at  A).  The  angle  from 

>  0  to  >  0  in  the  lower  coating  is  also  30®  but  is  90®-a  =  60®  in 
the  upper  coating. 

Measurements  were  made  at  points  with  integer  values  of  n^  from  1 
to  4,  with  a  bending  moment  causing  values  of  of  1  or  2  in  the  0.108 

in.  coating,  but  only  nj^  =  1  in  the  0.057  in.  coating  to  avoid  permanent 
straining, and  value  of  angle  a  from  0®  by  tens  to  90®.  At  first  a  point 
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of  selected  values  n^,  a  was  identified  on  the  lower  coating  by  cross¬ 
hairs  conveniently  fixed  in  a  wide  ring.  Figure  14a  shows  the  cross-hairs 
at  a  point  with  n^  =  2  (nj^  =  0)  in  plane  polarized  light  at  a  plate 
azimuth  making  its  neighborhood  brightest.  Fig.  14b  shows  the  same  point 
to  be  on  the  30®  isoclinic  (o  =  30®).  The  bending  moment  causing  nj^  =  1 
was  then  applied  and  the  fringe  order  n^  at  the  cross-hairs  was  estimated 
at  the  plate  azimuth  making  its  neighborhood  brightest  (Fig.  15a,  n^  =  1.9). 
The  angle  ^^2  “  “®u2  found  from  the  plate  azimuth  making  the  region  of 

the  cross-hairs  darkest  (Fig.  15b,  41^2  =  12®).  The  bending  was  then  in¬ 

creased  to  n,  =  2  and  the  new  values  of  n  and  4  _  were  measured  (Fig. 
b  u  u2 

16a;  =  2.8  and  Fig.  16b;  ^^2  =  After  all  measurements  on  the  lower 

coating  were  completed,  the  plate  was  turned  around  at  similar  measurements 
which  were  made  on  the  upper  coating. 

Test  Results 

The  results  are  surprisingly  good.  The  observed  total  fringe  order  n^ 

(cibove)  and  n^  (below)  and  their  "principal  directions"  are  shown  on  Mohr's 

circle  constructions  in  Figures  17  to  19.  As  shown  in  the  inset  to  Fig.  17 

the  s!J5)erposition  of  nj^  at  an  angle  a  to  n^  gives  the  total  fringe 

order  n^  =  NA’  for  the  coating  below.  For  other  angles  a  the  locus  of 

points  A'  will  be  a  circle  of  the  same  radius  n^  with  center  0  at  a 

distance  MO  =  n,  from  M  .  Likewise  for  other  values  of  n^  but  identical 
b  t 

n^^  the  total  fringe  points  (A')  will  lie  on  concentric  circles  with  0  as 
center  and  with  the  corresponding  n^  as  radii,  as  shown  in  Figures  17  to  19. 

The  magnitude  n^  is  found  as  the  distance  from  N  of  point  A'  determined 
by  the  angle  2a  at  0  . 


According  to  the  assumptions  made  in  the  present  paper  the  apparent 
principal  direction  should  be  the  one  at  the  free  surface ,  which  is  found 

(see  Fig.  5  for  the  stresses)  by  increasing  nj^  to  n^^H/a  =  AA":  the 

desired  direction  is  given  by  the  angle  ONA"  =  (fig*  17,  inset). 

However,  instead  of  drawing  a  new  set  of  concentric  circles  for  A"  ,  one 

has  only  to  take  the  same  distance  NF  =  A ’A"  =  n^^Ah/a  to  the  left  of 
N  for  all  values  n^  ,  and  determine  2(|»^2  the  angle  <  OFA'  .  This 
was  done  in  all  three  Figures  17  to  19.  The  total  fringe  order  n^2 
principal  direction  of  the  coating  above  is  found  in  the  same  way  by 

starting  with  an  angle  2a  -  180  at  0  (clockwise):  corresponding  pairs 
of  n^2  »  opposite  ends  of  the  same  diameter. 

In  Figures  17  to  19  the  experimental  points  have  been  plotted  by  draw¬ 
ing  a  line  through  F  at  angle  to  FN  and  taking  a  length  n^  from 

N  to  that  line.  The  experimental  points  have  been  joined  to  the  theoreti¬ 
cal  by  tail-like  segments ,  whose  magnitude  and  direction  indicate  the  corre¬ 
sponding  errors.  As  can  be  seen  the  results  are  extremely  good  except  when 
o  is  about  70°  or  80°  for  the  coating  below,  or  about  10°  to  20°  for  the 
coating  above,  as  was  expected  from  an  earlier  discussion  (Figures  3a,  b  and 
5).  Even  then  the  magnitudes  n^,  n^  are  almost  exact,  the  main  error  being 

in  the  directions  (b  ,  (t 

^u’  ^o 

A  computer  program  (both  program  and  tabulated  results  are  given  in 
the  Appendix)  was  set  up  for  the  calculation  of  n^  and  n^  corresponding 
to  the  total  stress  at  mid-thickness,  as  well  as  by  integration  (eq.  11),  and 


of  at  the  interface;  ^^2'  '^u2  free  surfaces;  and  (|)^^,  (j»^2 

at  mid- thickness .  The  directions  ,  (j)^2  ^t  the  free  surface  were  always 
closer  to  the  experimental  results.  Calculations  were  made  for  n^^  = 


1  and 
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=  2  with  the  0.108  in.  coatings  and  nj^  =  1  with  the  0.057  in.  coating. 

In  both  cases  the  values  of  n^  extended  from  1  to  4. 

The  same  program  gave  the  results  of  the  inverse  calculation  of  n^,  n^^, 
a  from  the  observed  values  n^,  n^,  <l>^  and  compared  them  with  their 

true  or  nominal  values.  Furthermore,  the  neglected  quantity  A  in  Eqs.  (17) 
and  (18)  and  the  resulting  percentage  error  in  6^  was  computed  at  each 
point.  The  error  was  throughout  either  exactly  zero  or  less  than  | 0.005 |  . 

The  nominal  total  rotation  through  the  coatings  below  and  above  and  the 

average  value  =  2(4»^2  “  *^ii^^”i  twice  the  rotation  over  the  coating 

thickness  to  the  retardation  for  a  stress  equal  to  the  one  at  mid-thickness 
were  also  computed.  These  values  are  plotted  in  Figures  3a,  b. 

The  results  of  the  inverse  calculation  (n^,  n,  ,  a  form  n  ,  n  ,  6  _, 

t’  b’  u’  o’  u2’ 

O^q)  are  shown  in  Figures  20  to  22,  where  they  are  compared  with  their  true 
values.  The  calculated  angle  a  is  taken  counterclockwise  from  the  hori¬ 
zontal  and  the  calculated  magnitude  of  n^^  is  marked  on  it  by  an  empty 
circle.  The  calculated  value  of  n^  is  taken  along  the  radius  correspond¬ 
ing  to  the  true  angle  a  and  marked  by  a  black  circle.  The  calculated 
points  are  joined  to  the  exact  points  by  tail-like  segments.  The  results 
are  remarkably  good  even  for  angles  a  of  eibout  70®  and  20®  where  respec¬ 
tively  6^2  substantial  errors.  No  systematic  pattern  of 

errors  is  apparent.  The  absolute  errors  appear  to  be  only  slightly  larger 
for  n^  =  nj^  =  1  (both  thicknesses)  than  for  higher  values,  but  the  rela¬ 
tive  errors  are,  of  course,  much  smaller  for  the  higher  values  of  n^  and 
njj  .  The  independence  of  absolute  error  from  the  magnitude  of  n^  and  n^^ 
may  be  the  result  of  a  constant  source  of  error,  such  as  the  small  accuracy 
of  measurement  (fringe  orders  measured  to  the  nearest  0.1  fringe,  angles  to 
the  nearest  1®),  which  was  purposely  chosen  so  as  to  simulate  practical 
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testing  conditions,  or  the  thin  cement  layer  (about  0.002  in.)  which  had 
no  "membrane"  stress.  However,  the  difference  between  observed  and  theo¬ 
retical  total  birefringence  is  largest  at  angles  a  around  70°  (error  in 
n^)  or  20°  (error  in  n^)>  Figures  17-19,  and  it  might  be  expected  that  the 
inverse  solution  for  n^^,  n^  might  also  be  less  accurate.  It  would  there¬ 
fore  be  useful  to  recognize  these  cases.  As  may  be  seen  in  Figures  17  to 
19  when  a  =  70°  or  a  =  20°  the  difference  In  -  n  |  is  large  and  the 
angle  2  1 0^^  ~  ®q2  I  “  ^  I  *^o2  "  ^u2  ^  zero .  For  zero  angle  it  is 

easily  found  that  a  =  0  or  as  90°  so  that  no  rotation  occurs  and  the 
results  should  be  of  the  highest  accuracy.  For  n^  “  n^,  o  should  be  close 
to  45°. 

It  is  true  that  in  the  present  tests  the  direction  of  the  bending  mo¬ 
ment  was  known  (major  diagonal)  and  did  not  have  to  be  calculated.  This 
calculation  may  easily  be  made  as  explained  at  the  end  of  the  previous 
paragraph  but  would  not  affect  the  calculated  values  of  n^^,  n^,  and  o  , 
as  the  only  possible  new  error  would  be  identical  in  both  directions 
and  Y^  of  M  and  N  and  not  in  their  relative  angle  a 

Conclusions 

The  tests  show  that  with  the  proposed  method  using  two  birefringent 
coatings,  the  principal  membrane  and  bending  stress  differences  can  be 
determined  with  an  accuracy  of  about  5%  to  10%  in  magnitude  and 
*2°  to  5°  in  direction.  This  compares  quite  favorably  with  the 
accuracy  (5%  and  2°)  expected  in  two-dimensional  photoelasticity. 
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APPENDIX 

Computer  Program  and  Tabulated  Results 

The  computer  program  is  given  in  pages  40-45  and  the  tabulated  results 
in  Tables  I  (0.108  in.  coatings)  and  II  (0.057  in.  coatings),  pages  46-52 
and  53-56  respectively. 

In  Tables  I  and  II  the  results  are  given  in  sets  of  3  lines,  each  set 
corresponding  to  a  single  point.  The  first  of  each  set  of  three  lines  (Row 
NOM)  gives  the  actual  or  nominal  values  of  the  fringe  order  n^  due  to  mem¬ 
brane  stress  only  (column  NT);  due  to  bending  only  (column  NB);  the 

angle  a  of  their  principal  directions  (column  ALFA);  the  total  expected 
fringe  orders  n^,  n^  (colvimns  NU,  NO)  below  and  above  and  surface  principal 
directions  41^,  <|)^  (columns  PU,  PO)  as  calculated  by  the  approximate  formu¬ 
las  (12)  and  (16). 

The  second  row  (EXP)  of  each  set  of  three  lines  gives  under  the  same 
colximns  the  experimental  values  n^,  n^,  «j)^,  (|)^  and  the  inversely  calculated 
values  of  n^,  nj^  and  a  according  to  formulas  (17,  with  A  =  0),  (13),  (7) 
and  (14). 

The  third  row  (INT)  gives  the  value  of  the  total  apparent  fringe  orders 
n^,  n^  (columns  NU,  NO)  as  found  from  the  integrated  form  (lla,b).  The  value 
of  the  neglected  quantity  A  computed  as  in  (18)  is  given  in  Row  INT  under 
column  PO;  and  the  percentage  error  in  6^  due  to  the  neglect  of  A  is 
given  under  column  PU.  As  can  be  seen  both  A  and  the  error  in  6^  are 
either  exactly  zero  (shown  by  0.0)  or  less  than  | 0.005 |  (shown  by  0.00). 

For  nj^=  n^  ,  a  =  0  ,  the  value  of  A  is  zero  but  the  percentage  error  in 
6^  is  indeterminate  (****),  The  last  but  one  column  DF2-1  gives  the  total 


rotation  of  principal  stress  in  degrees  from  free  surface  to  interface  in 
coating  below  (first  row)  and  above  (second  row).  The  last  column  RAV  E-3 
gives  the  value  x  10^  (radians  per  fringe)  of  the  average  rotation 

divided  by  half  the  total  retardation  in  the  coating  below  (first  row)  and 
above  (second  row). 
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COHPOTER  PROGRAM 


VEL  02  NOV.  66  OS/360  FORTRAN  H 

COMPILER  OPTIONS  -  NAME=  MAINt0PT=00,LINECNT»60«S0URCE»&C0tNQLlST»NQ0ECK»t0 


ISN 

0002 

DINENSION  AC2)9BC2)»CI2l»0C2)»Pi2l#Ef2l«FI2l«DF21i21t 
7  V(2l9UC2).Mr2l«JTAfil6}«XPI2)»XHINTt2l9SNC2)  9DF22i2) 

ISN 

0003 

READdtlOOOl  HlfH 

ISN 

0004 

1000 

FORMAT (4F10« 5) 

ISN 

0005 

H0=(HI+H)/2. 

ISN 

0006 

00  1423  I»lv6 

ISN 

0007 

1423 

JTABI 11=0 

ISN 

0008 

JTABC1}=1 

ISN 

0009 

JTABC5)=1 

ISN 

0010 

K=0 

ISN 

0011 

CALL  WRITE(HI,H) 

C  THE  NEXT  CARO  MUST  BE  CHANGED  IN  THE  FINAL  PROGRAM 

ISN 

0012 

JTAB(6)=0 

ISN 

0013 

1 

REA0(l9i0009EN0=1032}  SNU«SNOf RUtRO 

ISN 

0014 

THETAM=2.*CP0“RU) 

ISN 

0015 

THETAD=ABS ( THETAMI 

ISN 

0016 

THETA=  C  THETAM+3. 14159) /ISO.O 

ISN 

0017 

IF(91*0-THETA0)  30,20t20 

ISN 

0018 

20 

IF(THETAI>«89*0)  30,40,40 

ISN 

0019 

30 

EM==  ( ( SN0*^2-SNU**2 )  /  C  SN0**2+SNU**2 ) )  ♦« 

ISN 

0020 

SEM=$QRTC1.0“EM) 

ISN 

0021 

J^l 

ISN 

0022 

IF(EM-0.90)  13,13,11 

ISN 

0023 

13 

0EN=1.~(  l.-EM*IHQ/HI  ( COS  1  THETA l<^2 )  ) 

ISN 

0024 

IL=0 

ISN 

0025 

DELT2=I  C  l.+SEM^COSITHETA) 

ISN 

0026 

DELT»SQRT(0£LT2) 

ISN 

0027 

GO  TO  50 

ISN 

0028 

40 

DELT=1,0 

ISN 

0029 

J=-l 

ISN 

0030 

60  TO  50 

ISN 

0031 

11 

DELT*HI/HO 

ISN 

0032 

IL=1 

ISN 

0033 

50 

TEN=SQRT(C0.5*<SN0**24SNU^*2)l/ll*>l«DELT4HOJ/KTJ«^2l) 

ISN 

0034 

B  EN= ( T  EN^OELT  *H0 ) /HI 

ISN 

0035 

C0S2 A= ( SNU^*2-SN0**2 ) / 1 4* ♦TEN^BEN ) 

ISN 

0036 

OV=THETA 

ISN 

0037 

IF(CQS2A-1«I  49,18,18 

ISN 

0038 

18 

C0S2A=1,00 

ISN 

0039 

49 

IFC-C0S2A--1.)  97,98,98 

ISN 

0040 

98 

C0S2A=“1*00 

ISN 

0041 

97 

CONTINUE 

ISN 

0042 

ALFTR=ARC0SCC0S2A) 

ISN 

0043 

IF( J)  60,61,61 

ISN 

0044 

61 

ALFR=ALFTR 

ISN 

0045 

GO  TO  70 

ISN 

0046 

60 

IFCOV)  62,63,63 

ISN 

0047 

62 

ALFR=ALFTR 

ISN 

0048 

GO  TO  70 

ISN 

0049 

63 

AtFR=-ALFTR 

ISN 

0050 

70 

ALF^CALFR4'180,0)/I2«^3.  14159) 

ISN 

0051 

DH=HI-H 

ISN 

0052 

TZN=JTABI1) 

ISN 

0053 

RHO=  C JTAB( 6 ) ♦B. 141592)790. 

ISN 

0054 

SIN2A=SIN«RHQ) 

ISN 

0055 

C0S2A=C0SCRH0) 

ISN 

0056 

DELT={ JTAB(5)*HI )/(JTA8«l )*H0) 

ISN 

0057 

A(1)=DELT+C0S2A 

ISN 

0058 

A(2)=DELT-C0S2A 

-41- 


ISN 

0059 

B  Cl 1 » ( 0ELT«H )  /Hl^COSZ A 

ISN 

0060 

B  C  2)  DELT^HI  /HI--C0S2A 

ISN 

0061 

CC 1)*C0ELT^HQI/H1+C0S2A 

ISN 

0062 

C( 2)=C  DELT#H0>/HI-C0$2A 

ISN 

0063 

0(1)=C0S2A4>0ELT 

ISN 

0064 

D(2)=C0S2A-0ELT 

ISN 

0065 

E  C  l>*C0S2A+i0ELT^H)  Ml 

ISN 

0066 

E  ( 2 )  ^0S2A-  CDELT*H1/HI 

ISN 

0067 

DO  44  I»lt2 

C 

HERE  l*Uf  2=0 

ISN 

0068 

IF(ABSISIN2AH0»005)  89«d9«88 

ISN 

0069 

89 

XN1NT(1I=JTAB(1I<I'JTA&C5I 

ISN 

0070 

XN1NTI2)=JTAB(11-JTAB(5) 

ISN 

0071 

SNCI>=SQRTilT2N*A2>^tClil^^24SIN2A^^2n 

ISN 

0072 

GO  TO  87 

ISN 

0073 

88 

CONTINUE 

ISN 

0074 

V(I)=CACll^SQRTiAII)^«2^SIN2A«^2M/iBUHSQRTIBII}«^2^SIN2A4^#2ll 

ISN 

0075 

U(iy=ACI>«SQRTiACl|4^2^SlN2A«*2} 

ISN 

0076 

WC 1 I^BCl I^SQRTIBC 1 >«^2+SlN2A««2) 

ISN 

0077 

SN C I  >=SQRT  I C TZN«2 1  ♦! C ( I  l^Z^S I N2A#*2 1 1 

ISN 

0078 

V(l>=ABSCVtin 

ISN 

0079 

lFCABS(VCni-0«005)  33«33f32 

ISN 

0080 

33 

XNINTC  1)=^JTAB(  ll-JTABfS) 

ISN 

0081 

XNINT(2}=aTABtlKJTABC5) 

ISN 

0082 

GO  TO  87 

ISN 

0083 

32 

CONTINUE 

ISN 

0084 

XNINTCI)  =  (HI/(2.«0ELT4FDH)J«TZN*(UCI)-lillHISIN2A*«2l#AL0GI\ICIII  ) 

ISN 

0085 

87 

CONTINUE 

ISN 

0086 

F  (  U=ATAN2  ( S IN2 A «  E 1 1  n 

ISN 

0087 

P(I)=ATAN2(SIN2A«DI1)) 

ISN 

0088 

CALL  SlHZCFCDtRHOfREStll 

ISN 

0089 

F(I)==RES 

ISN 

0090 

CALL  SIHZCP(l>>RHO»RESf 11 

ISN 

0091 

PCD^RES 

ISN 

0092 

DF211II=PCI1-FIII 

ISN 

0093 

DF22 ( 1 1=  C  0F21 i I )«1S0. ) /3. 141592 

ISN 

0094 

PI  I )=( PC  I I«90.0l/3* 14159 

ISN 

0095 

FC I l=C  F C 1 }«90. 1/3. 14159 

ISN 

0096 

DF22CI>=DF22Cl>/2« 

ISN 

0097 

44 

CONTINUE 

ISN 

0098 

IFCSNU)  80«81t80 

ISN 

0099 

81 

RAVCl) =999999999*9 

ISN 

0100 

GO  TO  83 

ISN 

0101 

80 

CONTINUE 

ISN 

0102 

RAVC 11=1000. ♦0F21C1>/SNU 

ISN 

0103 

83 

IFCSNO)  d4»85984 

ISN 

0104 

85 

RAVC  2) =999999999. 9 

ISN 

0105 

GO  TO  86 

ISN 

0106 

84 

CONTINUE 

ISN 

0107 

RAVC21=1000.0»DF21C2)/SNQ 

ISN 

0108 

86 

CONTINUE 

ISN 

0109 

0H2=0H^*2 

ISN 

0110 

AA=1.-EM-“EM*CSINCTHETA)^*2)4C0H2/I8.0»HQ4H1)1^*2 

ISN 

0111 

IFCILI  300,300,301 

ISN 

0112 

301 

IFCAA) 302,303,303 

ISN 

0113 

303 

AC=2.*ISQRTCAA)-SEM) 

ISN 

0114 

GO  TO  132 

ISN 

0115 

302 

IFCAA^O.OOl)  130,305,305 

ISN 

0116 

305 

AA=0. 

ISN 

0117 

GO  TO  130 

ISM  0118  300  CONTINUE 

ISN  0119  IFCAAI  130tl31»131 

ISN  0120  131  A8^C 1.4SEM^C€SITH£TA1 

ISN  0121  AN==f2.>»CQ$ITHETAI|4t(SQRTlAAI-SEN) 

ISN  0122  IF(AB)  54,55t54 

ISN  0123  55  AC=-AA 

ISN  0124  GO  TO  56 

ISN  0125  54  AC=AN/AB 

ISN  0126  56  AA^AN 

ISN  0127  GO  TO  132 

ISN  0128  130  AC«999999*9 

C  THIS  NUMBER  IS  >THAN  THATAtLOWEO  BY  THE  FORMAT  AND  WILi.  PRINT  AS 
ISN  0129  132  CONTINUE 

ISN  0130  IF(K«14)  9lr92t92 

ISN  0131  92  WRITEC3f501) 

ISN  0132  K^^l 

ISN  0133  91  CONTINUE 

ISN  0134  IF( JTABClI-ll  220t41l9220 

ISN  0135  411  IFf  JTABC6l-*0)  22Qt2l0w220 

ISN  0136  210  WRITE(3tSQll 

ISN  0137  220  WRITE!  39502}  JTAB(  I }  9  JTABf  5)  9  JTA&I6J  tSNC  1)  t  SNUI  tPil )  t PI2)  9  0F22C  1 

?  9RAV(1) 

ISN  0136  230  WRITE! 39504}  TEN9BEN9ALF9 SNU9SN09PU9P09DF22!2I 9RAyi2l 

ISN  0139  240  WR1TE!39503}  XNINTC1)9XNINT(2}9AC9AA 

ISN  0140  501  FORMAT! •1*9 10X9*  NT  NB  ALFA  NO  NO  PU  P 

?  DF2-1  RAV  E“3*l 

ISN  0141  502  FORMAT! IHO9 10X9  «NOM* 9 1492X9 1492X9l6t3X9F4^29BX9F4»292X9F6«292X9 

?  F6.293X9F6.292X9F6.2i 

ISN  0142  503  FORMAT! IH  9iOX9»lNT  —  ~  —  • 9F4.293X9F4.293X9F6.2 

7  2X9F6.2} 

ISN  0143  504  F0RMAT!1H  9lOX9«EXP  •«  F4.292X9F4.292X9F6.292X9F4.293X9F4.292X9 

7  F6.292X9F6.293X9F6.292X9F6.2} 

1 SN  0 144  KTAB» J  T AB ! 6 } 

I SN  0 145  CALL  VINPT ! JTA89 KTAB } 

ISN  0146  K*K+1 

ISN  0147  GO  TO  1 

ISN  0148  1032  STOP 

JSN  0149  END 
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OS/360  FORTRAN  H 


COMPILER  OPTIONS  -  NAII£=  MAINf OPT^OOvLlNECNT^AQrSOURCEt BCO«NOLlSTtNODECKtLO 


ISN 

0002 

ISN 

0003 

ISN 

0004 

ISN 

0005 

ISN 

0006 

ISN 

0007 

ISN 

0008 

ISN 

0009 

600 

ISN 

0010 

601 

ISN 

0011 

602 

ISN 

0012 

603 

ISN 

0013 

ISN 

0014 

SUBROUTINE  NRITEIHI^MI 

WRlTEf3t6001 

URITE(3f601) 

URITE(3f602> 

A=2.6H 

B=HI-H 

NRITEC3«6031  AtB 

FORMAT i* I »»///////////////) 

FORMAT! iHatTOX^rSEPARAI ION  OF  MEMBRANE  AND  BENDING*//  llXt*STRES 
?  DIFFERENCES  IN  SHELLS  III TH*//  lAXr^TNQ  BtREFRINGENT  COATINGS* I 
FORMAT  ilM0» /////// /////yy//// 1 

FORMATlIHOt  1QX#*PLATE  THICKNESS  iINCHESI*«cF6^4/IIIX»«CQAIING  THI 
7KNESS  C  INCHEST  *:^6«4I 
RETURN 
END 
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VEL  02  NOV.  66  05/360  FORTRAN  H 

COMPILER  OPTIONS  -  NAME=  NAlN»OPT300,LINECNT960»SOIIRCE»BCO»NOtlST#NOOECK»JtO 


ISN 

0002 

SUBROUTINE  VINPTI JTABtKTABI 

ISN 

0003 

DIMENSION  JTA&f6l 

ISN 

0004 

IF(KTAB-40I  30,31,30 

ISN 

0005 

31 

IFCJTAB(l)-4}  32,33,33 

ISN 

0006 

33 

JTAB(6)^45 

ISN 

0007 

JTABfD^l 

ISN 

0008 

RETURN 

ISN 

0009 

32 

JTAB(1I=JTAB(1H1 

ISN 

0010 

RETURN 

ISN 

0011 

30 

IFIKTAB-45)  40,41,40 

ISN 

0012 

41 

lF(JTAB(ll-4}  42,43,43 

ISN 

0013 

42 

JTAB(1)=JTABMU1 

ISN 

0014 

JTAB(6)=45 

ISN 

0015 

RETURN 

ISN 

0016 

43 

JTABf6}=50 

ISN 

0017 

JTABClI^l 

ISN 

0018 

RETURN 

ISN 

0019 

40 

CONTINUE 

ISN 

0020 

IF(KTAB~90)  10,11,11 

ISN 

0021 

10 

CONTINUE 

ISN 

0022 

GO  TO  13 

ISN 

0023 

11 

IFIJTABC 11-41  13,21,21 

ISN 

0024 

21 

JTABC6)^0 

ISN 

0025 

JTABCll^l 

ISN 

0026 

JTAB(5)~JTAB(5H1 

ISN 

0027 

RETURN 

ISN 

0028 

13 

CONTINUE 

ISN 

0029 

JTAB(1I«:  JIABfl|4l 

ISN 

0030 

IFCJTABf 11-51  997,999,999 

ISN 

0031 

997 

RETURN 

ISN 

0032 

999 

JTABClI^l 

ISN 

0033 

12 

JTAB(6I=KTAB^10 

ISN 

0034 

RETURN 

ISN 

0035 

END 
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OS/360  FORTRAN  H 

COMPILER  OPTIONS  -  NAME=  MAIN>QPT«00,LlNECNT«60»SOURCE»BCO«NOLlSTtNODECK»LO 


ISN 

0002 

SUBROUTINE  SIHICSlt S2f S3f II 

ISN 

0003 

IFII-II  150.150tl51 

ISN 

0004 

150 

IF(S2)  152»153rl53 

ISN 

0005 

153 

IFCSl)  154tl55«T55 

iSN 

0006 

155 

S3«S1 

ISN 

0007 

RETURN 

ISN 

0008 

154 

S3=Sl+3. 14159 

ISN 

0009 

RETURN 

ISN 

0010 

152 

IFiSlI  156rl57,I57 

ISN 

0011 

157 

S3=:S1~3. 14159 

ISN 

0012 

RETURN 

ISN 

0013 

156 

S3*S1 

ISN 

0014 

RETURN 

ISN 

0015 

151 

IFCS2)  162tl63,163 

ISN 

0016 

163 

IF(Sl)  164tl659l65 

ISN 

0017 

165 

S3==Sl-3. 14159 

ISN 

0016 

RETURN 

ISN 

0019 

164 

S3=^S1 

ISN 

0020 

RETURN 

ISN 

0021 

162 

IFCSl)  166«1679l67 

ISN 

0022 

167 

S3=:S1 

ISN 

0023 

RETUlUC 

ISN 

0024 

166 

S3=:Sl+3.  14159 

ISN 

0025 

RETURN 

ISN 

0026 

END 
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TABLE  I 

TABULATED  RESULTS 
COATING  THICKNESS  0.1079  in. 
PLATE  THICKNESS  0.2470  in. 
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NT 

NB 

ALFA 

NU 

NO 

NQM 

1 

1 

0 

2.00 

0.00 

EXP 

1.00 

1.00 

0.0 

2.0U 

0.0 

INT 

“ 

— 

— — 

2.00 

0.0 

NOM 

2 

1 

0 

3.00 

1.00 

EXP 

2.00 

1.00 

0-0 

3.00 

1.00 

INT 

" 

— 

— — 

3.00 

1.00 

NOM 

3 

1 

0 

4.00 

2-00 

EXP 

3.00 

1.00 

0.0 

4.00 

2-00 

INT 

— 

4.00 

2.00 

NOM 

4 

1 

0 

5.00 

3.00 

EXP 

4.00 

1.00 

0.0 

5-00 

3-00 

INT 

— — 

— 

— — 

5.00 

3.00 

NOM 

1 

1 

10 

1.97 

0.35 

EXP 

1.19 

0.94 

8.47 

2.10 

0.40 

INT 

— — 

1.97 

0.39 

NOM 

2  . 

1 

10 

2.96 

1.11 

EXP 

2.04 

0.97 

5.53 

3.00 

1.10 

INT 

— 

— 

2.96 

1.12 

NOM 

3 

1 

10 

3.95 

2-09 

EXP 

3.04 

0.98 

6.73 

4.  00 

2-10 

INT 

—  — 

— — 

3.95 

2.09 

NOM 

4 

1 

10 

4.95 

3.08 

EXP 

3.98 

1.06 

9.64 

5.00 

3-00 

INT 

“ 

— 

— —  . 

4.95 

3.08 

NOM 

1 

1 

20 

1.88 

0.68 

EXP 

1.18 

0.92 

18.19 

2.00 

0.70 

INT 

— ■ 

— 

1.88 

0.70 

NOM 

2 

1 

20 

2.84 

1.39 

EXP 

2.06 

0.97 

18.18 

2.90 

1.40 

INT 

— — 

— 

— — 

2.84 

1.40 

NDM 

3 

1 

20 

3.82 

2.32 

EXP 

3.02 

1.06 

19.49 

3.90 

2.30 

INT 

— 

3.82 

2.33 

NOM 

4 

1 

20 

4.81 

3.30 

EXP 

3.99 

1.08 

18.57 

4.90 

3.20 

INT 

—  — 

— — 

4.81 

3.30 

NOM 

1 

1 

30 

1.73 

1.00 

EXP 

1.14 

1.01 

27.62 

1.90 

1.00 

INT 

“ 

— 

1.73 

1.01 

NOM 

2 

1 

30 

2-65 

1.73 

EXP 

2.18 

1.04 

27.66 

2.90 

1.60 

INT 

— 

— 

— 

2.65 

1.74 

PU 

PO 

DF2-1 

RAV  E-3 

0.0 

0.0 

0.00 

-90.00 

0.0 

0.0 

90.00 

0.0 

«****«: 

0.0 

0.0 

0.0 

-90.00 

-90.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-90.00 

-90.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-90.00 

-90.00 

0.0 

0.0 

0.0 

0.0 

0.0 

4.33 

5.00 

-0.00 

-21.60 

-38.00 

-0.00 

-1.57 

41.14 

-26.12 

6.06 

7.00 

0.00 

-65.03 

-67.00 

0.00 

-1.37 

9.95 

-15.99 

315.80 

6-99 

8.00 

0.00 

-72.95 

-73.00 

0.00 

-1.15 

4.16 

-10.03 

69.10 

7-56 

7.00 
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FIG.  I  DOUBLY  COATED  SHELL  IN  SIMPLE  TENSION  AND  IN  PURE 
BENDING. 
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FIG.  2  PRINCIPAL  DIRECTIONS  OF  MEMBRANE,  BENDING,  AND 
TOTAL  STRESS  IN  THE  LOWER  COATING,  SEEN  FROM 
BELOW  Z  >  0. 
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FIG.  3a  AVERAGE  RATE  OF  ROTATION  (Ravg)  'N  LOWER 
COATING  vs.  THE  ANGLE  a  BETWEEN  MEMBRANE 
AND  BENDING  OF  VARIOUS  INTENSITIES,  COATING 
THICKNESS  0.  lOBin. 


FIG. 3b  AVERAGE  RATE  OF  ROTATION  {  Ravg)  IN  LOWER 
COATING  vs.  THE  ANGLE  a  BETWEEN  MEMBRANE 
AND  BENDING  OF  VARIOUS  INTENSITIES.  COATING 
THICKNESS  0.057in. 


FIG.  4.  MOHR  DIAGRAM  OF  THE  SUPERPOSITION  OF 
MEMBRANE  STRESS  CTf  AND  BENDING  STRESS 
DIFFERENCE  crt,z  AT  DISTANCE  z>0  (COATING 
BELOW)  AND  -z  (ABOVE). OBSERVATION  IS  ALWAYS 
ALONG  NORMAL  TO  SHELL  FROM  z>0  TOWARD 
z<  0. 


COATING  ABOVE 


FIG.  5  STRESS  VARIATION  THROUGH  COATING 

THICKNESS  ROTATION  (^uZ'^^UI  (BELOW) 
AND  (i>o2-4>a\  (ABOVE).  Rmox  OCCURS  AT 
POINT  C.!  OF  LOWER  COATING 


M/M,  =  njj(  FRINGES)  n^(  FRINGES) 

G.  8.  VALUES  OF  ,  n  ^  PRODUCING  ONE  FRINGE  AT  Q  «  0»  I5®»30*.  45 
90®  AND  ONE -HALF  FRINGE  AT  0=60®  AND  75®. 


FIG.  9  SQUARE  PLATE  WITH  BIREFRINGENT  COATINGS 
ON  BOTH  FACES  SUBJECTED  TO  ANTICLASTIC 
BENDING.  DIAL  AT  LEFT  BELONGS  TO  DYNAMOMETER 
MEASURING  CORNER  LOADS.  DIAL  AT  CENTER 
HANGS  FROM  DEFLECTOMETER  SPANNING  FRONT- 
TO-BACK  DIAGONAL. BOTTOM  RIGHT  KNOB  IS  FOR 
LOAD  APPLICATION. 


FIG.  10  SCHEMATIC  DIAGRAM  OF  APPARATUS 


'  '’’W 


GENERAL  VIEW  OF  POLARISCOPE  AND  PLATE  IN  ANTICLASTIC 
BENDING  OBSERVED  AT  NORMAL  INCIDENCE. 


FIG.  12 

FRINGE  PATTERNS 
IN  LOWER  COATING 


a.  MEMBRANE  STRESS 
ONLY.  PRINCIPAL 
DIRECTIONS  ARE 

circumferential: 

nt;nb=o 


b.  BENDING  CAUSING 
ONE  FRINGE  SUPER¬ 
IMPOSED  ON  MEMB¬ 
RANE  STRESS. 
PRINCIPAL  BENDING 
DIRECTION  VERTICAL. 


C.  BENDING  CAUSING 
TWO  FRINGES  SUPER¬ 
IMPOSED  ON  MEMB¬ 
RANE  STRESS. 
PRINCIPAL  BENDING 
DIRECTION  VERTICAL. 

n^; 2 


FIGJ3.  POSITION  OF  CORRESPONDING  POINTS  P',  P"  IN  COATINGS 
BELOW  AND  ABOVE  AND  ANGLES  BETWEEN  POSITIVE 
MEMBRANE  AND  BENDING  STRESSES. 


FIG.  14  DETAIL  OF  PLATE  WITH  MEMBRANE  STRESS 
ONLY.  CROSS-HAIRS  AT  POINT  WITH 
n^=2,nj,=  0  AND  PRINCIPAL  DIRECTION  AT 
30®  TO  VERTICAL  DIAGONAL. 

COMPARE  WITH  FIG.  12  a. 

a.  PLANE  POL.  LIGHT  AT  45“  TO  PRINCIPAL  STRESS 

b.  PLANE  POL.  LIGHT  PARALLEL  TO  PRINCIPAL  STRESS 


FIG.  15  THE  SAME  POINT  AS  IN  FIG.  14  WITH  A 
SUPERIMPOSED  BENDING  CAUSING  ONE 
FRINGE:  n|  =  2,  =  I. 

COMPARE  WITH  FIG.  12  b.  Hy  =  1.9. 

a.  PLANE  POL.  LIGHT  AT  45®  TO  AZIMUTH  OF  MINIMUM 
INTENSITY 

b.  PLANE  POL.  LIGHT  AT  AZIMUTH  OF  MINIMUM  INTENSITY. 


FIG.  16  THE  SAME  POINT  AS  IN  FIG.  14  WITH  A 
SUPERIMPOSED  BENDING  CAUSING  TWO 
FRINGES*  nt  =  2,  n|3  =  2. 

COMPARE  WITH  FIG.  12b.  Py  =  2.8. 

a.  PLANE  POL.  LIGHT  AT  45*  TO  AZIMUTH  OF  MINIMUM 
INTENSITY. 

b.  PLANE  POL.  LIGHT  AT  AZIMUTH  OF  MINIMUM  INTENSITY. 


FIG. 17.  OBSERVED  ^ND  no,<^02  PLOTTED  ON  MOHR'S 

DIAGRAM  FOR  COMPARISON  WITH  THEORETICAL  VALUES. 
COATING  THICKNESS  Ah  «  0.108  IN.  0^,=  I  •,  «  I  TO  4  . 


FfG.  18.  OBSERVED  AND  no.</>o2  PLOTTED  ON 

MOHR'S  DIAGRAM  FOR  COMPARISON  WITH 
THEORETICAL  VALUES.  COATING  THICKNESS 
Ah  «  0.108  IN.  nj,*2;  *1  TO  4. 


FIG. 19.  OBSERVED  ny  AND  no,<^o2  PLOTTED  ON 

MOHR'S  DIAGRAM  FOR  COMPARISON  WITH 
THEORETICAL  VALUES.  COATING  THICKNESS 
Ah  ?  0.057  IN.  njj*l  -,  *  I  TO  4. 


FIG.  20  SOLUTION  FOR  THE  STRESS  -  DIFFERENCES. 

EXPERIMENTAL  RESULTS  SHOWN  BY  DOTS 
JOINED  TO  TRUE  VALUE  FOR  COMPARISON 
COATING  Ah  «  0.108  IN.; 


FIG.  21  SOLUTION  FOR  THE  STRESS  -  DIFFERENCES. 

EXPERIMENTAL  RESULTS  SHOWN  BY  DOTS 
JOINED  TO  TRUE  VALUE  FOR  COMPARISON. 
COATING  Ah  *0.108  ini 


n^*  I  0^*3 


FIG.  22  SOLUTION  FOR  THE  STRESS- DIFFERENCES 
EXPERIMENTAL  RESULTS  SHOWN  BY  DOTS 
JOINED  TO  TRUE  VALUE  FOR  COMPARISON 
COATING  Ah -0.057  in;  n,^  ?  I . 


